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ABSTRACT 


This  report  deals  with  pre-  and  post-test  work  on  Air  Force 
Structures  Program,  Project  5.3,  Tests  on  the  Loading  of  Horizontal 
Cylindrical  Shapes.  The  general  objective  of  this  test  was  to 
Increase  the  knowledge  of  blast  loadings  on  structures  of  cylindrical 
shape.  The  specific  objectives  were  to  determine  the  effect  of 
(l)  shock  strength,  (2)  size  (i.e.,  scaling),  (5)  ground  proximity, 
and  {k)  three-dimensionality. 

Five  horizontal  cylinders  of  various  size  and  distance  above 
ground  were  situated  in  the  Mach  region  of  both  Shots  9 10. 

Instrumentation  consisted  of  pressure  gages  and  strain  gages,  the 
latter  representing  the  output  of  a net  force  measurement  system. 

Most  of  the  test  objectives  could  not  be  realized.  The  majority 
of  the  pressure  data  exhibited  a considerable  and  generally  random 
behavior  which  obscured  the  effects  being  studied.  In  addition,  the 
strain  data  did  not  yield  any  quantitative  results;  however,  it  appears 
to  be  mainly  due  to  an  inability  to  interpret  the  data  properly. 

While  the  methods  of  data  reduction  employed  have  been  rather  unsuc- 
cessful, they  by  no  means  constitute  an  exhaustive  study  of  this 
problem.  Accordingly,  it  was  not  possible  to  reach  a definite  con- 
clusion as  to  the  feasibility  of  net  force  measurement  systems  whose 
output  represent  the  response  of  a complicated  dynamic  system  to  an 
input  which  is  itself  the  desired  information. 

The  pressure-time  distribution  was  obtained  on  a cylinder  one 
radius  above  the  ground.  The  results  indicate  the  build  up  and  clear- 
ing times  to  be  essentially  the  same  as  for  a cylinder  in  free  space, 
except  for  a somewhat  longer  clearing  time  on  the  front  face. 

The  analysis  of  pressure  gage  data  obtained  from  a vertical 
cylindrical  storage  tank,  included  as  part  of  Air  Force  Project  5.26.1. 
Test  of  the  Effects  on  POL  Installations  (WT-736),  is  presented  in 
this  report.  These  results  Indicate  that  the  finite  size  of  the  tank 
did  not  affect  the  front  face  loading  appreciably;  the  loading  on  the 
back  side  built  up  somewhat  more  rapidly  than  predicted  for  a two- 
dimensional  cylinder. 


iiCKW  neaTtnefCP 


FOREWORD 


This  report  is  one  of  the  reports  presenting  the  results  of  the 
78  projects  participating  in  the  Military  Effects  Tests  Program  of 
Operation  UPSHOT -KNOTHOLE,  which  included  11  test  detonations.  For 
readers  interested  in  other  pertinent  test  information,  reference  is 
made  to  WT-782,  Summary  Report  of  the  Technical  Director.  Military 
Effects  Program.  This  summary  report  includes  the  following  infonna- 
tion  of  possible  general  interest. 

a.  An  over-all  description  of  each  detonation,  including 
yield,  height  of  burst,  ground  zero  location,  time  of 
detonation,  ambient  atmospheric  conditions  at  detona- 
tion, etc.,  for  the  11  shots. 

b.  Compilation  and  correlation  of  all  project  results  on 
the  basic  measurements  of  blast  and  shock,  thermal 
radiation,  and  nuclear  radiation. 

c.  Compilation  and  correlation  of  the  various  project 
results  on  weapons  effects. 

d.  A summary  of  each  project,  including  objectives  and 
results. 

e.  A complete  listing  of  all  reports  covering  the 
Military  Effects  Tests  Prograun. 
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PREFACE 


In  a letter  dated  12  March  1952,  the  Air  Materiel  Connnand  was 
requested  by  Air  Research  and  Development  Command  to  submit  for  testing 
in  Operation  UPSHOT -KNOTHOLE  existing  requirements  for  a structures 
progreun  which  would  be  based  on  the  needs  of  the  Air  Force  for  target 
analysis  and  Indirect  bomb  damage  assessment  Information.  Within  the 
Air  Materiel  Command,  the  responsibility  for  designing  and  executing 
such  a program  was  delegated  to  the  Special  Studies  Office,  Engineering 
Branch  of  the  Installations  Division.  The  requirements  which  were 
submitted  and  approved  became  part  of  Program  3 of  the  operation  and 
were  designated  as  Projects  ^.1,  5*3»  3*5»  3*6  and  5.26.1. 

Mr.  B.  J.  O'Brien  of  the  Special  Studies  Office  was  appointed  Project 
Officer  and,  as  such,  coordinated  and  successfully  directed  the  plan- 
ning and  operational  phases  of  five  of  the  six  projects.  Due  to  the 
similarity  in  test  objectives  involving  railroad  equipment,  the  projects 
proposed  by  the  Transportation  Corps,  U.  S.  Army,  and  the  U.  S.  Air 
Force  were  combined  into  Project  5«6  with  Lt.  Colonel  Donald  G.  Dow, 

TC,  USA,  as  Project  Officer  and  Mr.  O'Brien  as  Assistant  Project  Officer. 

Armour  Research  Foundation  (ARF)  of  the  Illinois  Institute  of 
Technology  was  awarded  a contract  to  assist  the  Special  Studies  Office 
in  planning  and  designing  the  experiments,  and  in  analysis  and  report- 
ing of  test  results.  During  the  period  of  planning,  close  liaison  was 
maintained  with  other  interested  Air  Force  agencies,  particularly  the 
Physical  Vulnerability  Division,  Directorate  of  Intelligence,  Head- 
quarters, USAF.  Many  valuable  suggestions  were  contributed  by  Colonel 
John  Weltman,  USAF,  Lt.  Colonel  John  Ault,  USAF,  Messrs.  R.  G.  Grassy 
and  S.  White,  Dr.  F.  Genevese  and  others  of  that  Division,  and  by 
Mr.  Louis  A.  Nees,  Chief  of  the  Engineering  Branch,  Installations 
Division,  AMC. 

Personnel  of  the  Special  Studies  Section  who  were  intimately  con- 
nected with  the  program  were  Mr.  Eric  H.  Wang,  Chief,  Special  Studies 
Office,  who  was  the  technical  and  scientific  monitor  for  the  Air  Force 
Program,  Mr.  Arthur  Stansel,  and  Mrs.  Malsie  G.  Ridgeway,  secretary  to 
Mr.  Wang.  Other  members  of  the  Office  who  were  associated  with  the 
program  were  Messrs.  R.  R.  Blrukoff,  P.  A.  Cooley,  J.  C.  Noble,  and 
Lts.  T.  M.  Murray  and  G.  A.  Rockwell,  USAF. 


Most  of  the  Introduction  section  of  this  report  was  tak. n from 
the  preface  of  the  Preliminary  Report,  Operation  UPSHOT -KKOTilOLE, 
Project  5.5,  authored  by  Eric  H.  Wang  and  Bernard  J.  O'Brier. 

The  responsibility  within  the  Air  Force  for  execution  of  the  six 
projects  was  transferred  from  the  Special  Studies  Office,  Installations 
Division,  Air  Materiel  Command  to  Blast  Effects  Research,  Mechanics 
Branch,  Aeronautical  Research  Laboratory,  Wright  Air  Development  Center, 
on  15  November  195*<’* 
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1.1  PURPOSE  OF  AIR  FORCE  TEST  PROGRAMS 

The  series  of  tests  conducted  by  the  Air  Force  in  Operation 
UPSHOT -KNOTHOLE  is  part  of  a continuing  Air  Force  progreun  designated 
as  "Oetennlnatlon  of  Blast  Effects  on  Buildings  and  Structures."  The 
United  States  Air  Force  is  mainly  Interested  in  the  offensive  aspects 
of  such  research. 

The  UPSHOT -KNOTHOLE  projects  si>onBored  by  the  Air  Force  and  their 
specific  objectives  cannot  be  fully  understood  without  some  knowledge 
of  the  general  objectives  of  the  over-all  program.  The  research  results 
emanating  from  these  studies  and  experiments  conducted  by  the  Air  Force 
are  used  by  a number  of  government  agencies  to  improve  their  own  systems 
of  determining  blast  effects,  or  to  further  their  own  research. 

One  of  these  agencies  is  the  Directorate  of  Intelligence,  Head- 
quarters, USAF,  %rtiich  feeds  results  as  they  are  obtained  into  its  own 
system  of  vulnerability  classes,  thereby  making  it  possible  to  analyze 
prospective  enemy  targets  with  greater  accuracy,  and  to  recommend  the 
desired  ground  zero.  Another  principal  user  of  the  research  results  is 
the  Strategic  Air  Ccsnmand,  which  applies  them  toward  Improvement  of  an 
existing  Indirect  bomb  damage  assessment  system.  The  purpose  of  this 
system  is  to  make  it  possible  to  dispense  with  the  usual  reconnaissance 
after  a strike,  using  instead,  information  on  the  actual  ground  zero, 
hel^t  of  burst,  and  yield  of  the  weapon  which  is  brought  back  to  the 
operational  base  by  the  strike  aircraft  to  determine  the  damage  inflicted. 

The  task  of  determining  the  effect  of  blast  on  various  types  of 
building  structures  and  tactical  equipment  is  a rather  formidable  one. 
However,  its  difficulty  is  somewhat  relieved  by  the  fact  that,  for  the 
offensive  purposes  in  which  the  Air  Force  is  interested,  it  is  not 
necessary  to  determine  the  effect  of  transient  loads  on  these  items  with 
the  same  accuracy  as  would  normally  be  employed  for  static  design  pur- 
poses. In  fact,  even  if  it  were  possible  to  solve  the  dynamic  problems 
satisfactorily.  Intelligence  information  would  be  far  too  sketchy  to 
furnish  the  information  necessary  to  Justify  the  use  of  an  accurate 
analysis  for  items  located  in  prospective  enemy  countries.  From  the 
experience  that  is  so  far  available  it  is  expected  that  it  will  be 
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possible  within  the  foreseeable  future  to  determine  blast  daimage  within 
broad'  limits  with  sufficient  accuracy  for  planning  as  well  as  for  oper- 
ational purposes. 

In  view  of  the  complex  phenomena  attending  shock  waves  enianating 
from  various  types  of  atomic  blasts  and  the  uncertainties  inherent  in 
determining  significant  parameters,  an  investigator  s first  idea  would 
be  to  obtain  solutions  through  a long  series  of  very  elaborate  and 
properly  designed  full-scale  tests.  However,  neither  funds  nor  time 
will  allow  such  an  approach.  It  has  therefore  been  the  objective  of 
the  agencies  involved  to  obtain  sufficiently  accurate  results  by  Judi- 
cious use  of  theoretical  analyses,  laboratory  tests,  high-explosive 
field  tests,  and  a small  number  of  full-scale  atomic  tests. 

Three  of  these  research  projects  have  involved  full-scale  atomic 
testing.  The  first  was  (3EE21H0USQ,  the  second  was  JANCILB  (the 
first  imderground  burst  of  an  atoiaic  weapon  to  i^ilch  an  Air  Force 
structures  prograa  was  subjected)  and  the  third  the  UP3I0r-KN0TH0I£ 
program. 

From  previous  analysis,  laboratory  tests,  and  full-scale  tests 
(the  latter  especially  as  conducted  in  GREENHOUSE),  methods  of  damage 
prediction  have  been  developed  by  Armour  Research  Foundation  (ARF)  and 
others.  These  prediction  methods  have  attempted  to  describe  the  char- 
acter of  the  blast  loads  acting  on  a variety  of  items.  Response  compu- 
tations based  on  the  predicted  loadings  permit,  in  turn,  an  estimate  of 
physical  damage.  However,  the  relation  between  the  deflection  or  move- 
ment of  a body  and  significant  military  damage  has  never  been  clearly 
established  except  for  extreme  cases,  e.g.,  total  destruction  or  no 
destruction.  Another  aim  of  these  tests  is,  therefore,  to  establish 
the  relationship  between  deflection  and  functional  damage.  A full-scale 
test  also  affords  an  excellent  opportunity  to  determine  scaling  check 
points  for  laboratory  tests. 

In  addition  to  the  scientific  aspects  of  the  tests,  most  of  the 
results  of  the  Air  Force  projects  can  be  used  by  other  government  agen- 
cies such  as  the  Directorate  of  Intelligence  to  furnish  "rough  and 
ready"  experimental  answers  to  the  behavior  of  various  kinds  of  struc- 
tures under  blast.  In  many  cases  there  is  a statistically  significant 
number  of  items  involved  which,  added  to  previous  experimental  data 
such  as  those  gathered  at  Hiroshima  and  Nagasaki,  will  help  round  out 
the  present  vulnerability  picture.  In  other  cases,  mathematical  analysis 
may  have  to  rely  on  ad  hoc  information  to  furnish  parameters  which  cannot 
be  obtained  in  any  other  way. 

The  foregoing  remarks  are  designed  to  furnish  the  background  neces- 
sary for  a full  understanding  of  the  objectives  of  this  and  other  of  the 
Air  Force  projects.  The  full  significance  and  value  of  the  results  of 
each  test  will  be  realized  only  when  they  are  correlated  with  results 
of  past,  current,  and  future  analyses,  laboratory  tests,  high-explosive 
field  tests,  and  full-scale  atomic  investigations. 

1.2  SPECIFIC  OBJECTIVES 

Knowledge  of  the  blast  loading  on  circular  cylindrical  bodies  is 
required  in  connection  with  the  prediction  of  damage  to  certain  Items 
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of  target  interest.  Fuel  storage  tanks,  tank  cars,  chimneys,  and  pen- 
stocks are  but  a few  items  which  are  idealized  as  cylinders  for  the 
prediction  of  blast  loading. 

Existing  load-prediction  theories  are  based,  for  the  most  part, 
on  small-scale  model  studies  conducted  in  a shock  tube.  Thus,  a great 
deal  of  interest  concerns  the  validity  of  scaling  shock  tube  data  to 
full-scale  field  conditions.  Cylinders  of  target  interest  are  often 
supported  at  varying  distances  above  the  ground  or  other  reflecting 
surfaces.  Shock  tube  studies  of  this  configuration,  (Air  Force  Plan- 
ning ProgrEun  Final  Pretest  Report,  Part  I,  Oil  Storage  and  Equipment 
Tests )l/  have  indicated  an  increase  in  total  force  due  to  multiple 
reflections  between  the  undersurface  of  the  cylinder  and  the  ground. 

It  is  therefore  desirable  to  establish  whether  or  not  this  phenomenon 
occurs  under  field  conditions,  and  the  extent  to  which  it  may  be  pre- 
dicted. 

Since  most  shock  tube  studies  have  to  date  been  conducted  on 
two-dimensional  cylinders,  it  is  also  of  interest  to  compare  these 
results  with  the  loading  on  finite  (i.e.,  three-dimensional  cylinders). 

The  puiTpose  of  this  test  is  to  obtain  information  on  the  above 
items  by  determining  the  blast  loads  on  nondeformable  right  circular 
cylinders  supported  horizontally  at  various  distances  above  the  ground. 
The  specific  test  objectives  are  to  determine  the  influence  of  the 
following  parameters  on  the  blast  loading  of  cylinders: 

1 . Shock  strength 

2.  Size 

5 . Ground  proximity 

h.  Three-dimensionality 

The  aim  of  this  test  is  thus  to  obtain  information  which  can 
provide  the  proper  scale  relationships  for  continuing  shock  tube  and 
high-explosive  studies  on  the  blast  loading  of  circular  cylinders. 

This  information  can  be  ultimately  incorporated  in  suitable  damage 
prediction  schemes  for  items  of  this  type. 

1.3  RESPONSIBILITIES 

Armour  Research  Foundation  (ARF)  was  retained  by  the  Air  Materiel 
Command  (AMC)  of  the  United  States  Air  Force  to  carry  out  the  follow- 
ing work: 

A.  Consultation  on  the  selection  of  the  test  items. 

B.  Design  of  the  test  items. 

C.  Specification  of  instrumentation  requirements. 

D.  Location  of  the  structures  at  the  test  site. 

E.  Supervision  of  construction  of  the  tesft  items. 

"y  Complete  reference  concerning  publications  mentioned  in  parentheses 
may  be  found  in  the  Bibliography  at  the  end  of  this  report. 
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F.  Theoretical  and  experimental  analyses  concerning  pretest 
predictions  of  blast  loading  and  response  of  the  test  items 
where  required. 

G.  Analysis  of  the  test  results. 

H.  Submission  of  reports  accounting  for  the  Foundation's 
activities  pursuant  to  the  objectives  of  the  program. 

Detailed  statements  of  the  duties  and  obligations  of  the  contracting 
parties  can  be  found  in  the  Statement  of  Work  section  of  Air  Force 
Contract  AF55 (058) -30029. 

Preparation  of  the  construction  drawings  for  most  of  the  test 
items  was  subcontracted  by  the  Foundation  to  the  firm  of  Holabird  and 
I Root  and  Burgee.  A member  of  this  organization  supervised  the  actual 

r construction  under  the  general  direction  of  ARF.  As-built  drawings  of 

(all  the  test  items  were  prepared  by  the  Silas  Mason  Company. 

All  instrumentation  was  installed  and  operated  by  the  Ballistic 
Research  Laboratories  (BRL)  under  UPSHOT-KNOTHOLE  Project  5«28.1 
(Structures  Instrumentation,  WT-738).  BRL  was  also  responsible  for 
the  reduction  and  presentation  of  the  instrument  records. 


I 


16 


CHAPTER  2 


GENERAL  DESCRIPTION  OF  TEST 


2.1  TEST  ITEMS 

Five  horizontal  cylinders,  closed  at  both  ends  and  supported  at 
varying  distances  above  the  ground,  were  located  at  two  stations  in 
both  Shots  9 and  10.  Four  of  the  cylinders  were  identical  in  size  and 
shape;  the  fifth  was  a quarter-scale  model  of  the  other  four.  A sched- 
ule of  the  cylinder  dimensions  is  shown  in  Fig.  2.1;  as-built  construc- 
tion drawings  are  shown  in  Fig.  2.2. 

The  basic  cylinder  was  6 ft  in  dieuneter  and  20  ft  long.  Each 
cylinder  was  oriented  so  that  the  shock  wave  would  be  normal  to  its 
longitudinal  axis.  All  five  were  fabricated  from  l/2-in.  steel  plate 
with  reinforced  end  sections.  All  seams  were  butt-welded  so  that 
surfaces  of  interest  would  be  smooth.  The  cylinders  were  supported  on 
specially  designed  sensor  bars  welded  to  the  end  plates.  The  sensors, 
in  turn,  were  bolted  to  piers  of  reinforced  concrete  construction  set 
in  the  ground  (see  Fig.  2.2).!^  The  design  was  such  that  no  permanent 
deformation  or  gross  motion  of  the  cylinder  was  expected  at  the  intended 
overpressures.  Photographs  of  the  quarter-scale  cylinder  (5.3b)  and  an 
end  view  of  a typical  6-ft  diameter  cylinder  are  shown  in  Figs.  2.5  a^d 
2,4.  A summary  of  the  cylinder  properties  is  given  in  Table  2.1. 

An  analysis  of  the  pressure  gage  data  obtained  from  a vertical 
cylindrical  storage  tank,  included  as  part  of  UPSHOT-KNOTHOLE  Project 

5.26.1  (Tests  of  the  Effects  on  POL  Installations,  VIT-736),  is  presented 
in  this  report . The  tank  (3.26bd)  was  of  l/4-ln.  welded  steel  plate 
construction,  15  ft  in  diameter  and  10  ft  high.  The  roof  had  a slope 

of  approximately  10  degrees  and  consisted  of  a l/l6-in. - sheet  steel 
cover  on  twelve  l-by-4-in.  wooden  rafters.  The  rafters  were  toe-nailed 
to  a 4-by-4-in.  wooden  centerpost  and  were  strengthened  by  bridging. 

The  tank  was  filled  with  water  to  a level  of  7 ft.  A photograph  of 
the  tank  is  shown  in  Fig.  2.5. 
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T/iBLE  2.1  - Properties  of  Test  Cylinders 


Cylinder 

Property 

Units 

3 • 3^ f Of  d y 6 

3.3b 

Cross-sectional  area 

in.2 

U2 

27.5 

Mass  per  Unit  Length 
(Neglecting  End  Plates) 

lb  sec.2/in. 

0.0824 

0.0202 

Total  ISass  (including 

End  Plates) 

lb  sec. Via. 

28.5 

1.71 

Moment  of  Inertia  of 
Cross  Section  About 
Diameter 

in.‘^ 

1 

71.8  X ICP 

1.05  X ICP 

Stiffness  of  Sensor 
Supports* 

Ib/in. 

2.5  X lof 

5.5  X 10^ 

1.6  X 10^ 
2.7x10^ 

* "Stiffness"  refers  to  the  load  per  imit  deflection  of  the  sensor  bar 
assesibly  in  the  plane  of  the  sensor  bar  frame.  Due  to  the  geoiiietry 
involved,  the  stiffness  of  the  frame  is  Indepeisdent  of  the  orients^ 
tion  of  the  load  and,  hence,  is  given  by  the  load  per  unit  deflection 
of  a single  sensor  bar  loaded  axially  (stiffness  due  to  bending  of 
the  bars  is  neglected  throu^out).  The  larger  value  of  stiffness  is 
based  on  the  unsupported  sensor  bar  length;  the  lower  value  on  the 
total  length  between  center  of  the  cylinder  and  base  plate.  (The 
actual  value  of  stiffness  is  intermediate  to  these  two  extremes.) 

See  Fig.  2.2  for  dimensions  of  the  cylinders  and  sensor  frames. 

2.2  INSTRUMEMTATIQH 
2.2.1  General 

It  was  desired  to  measure  net  forces  on  the  cylinder  as  well  as 
local  forces  on  certain  portions  of  the  cylinder  (e.g.,  on  the  surface 
facing  the  ground).  Instrumentation  cozisisted  of  both  pressure  gage 
and  strain  gage  measurements.  Six  air-pressure  gages  were  mounted  on 
the  circumference  of  each  cylinder  at  mid-length.  The  gage  locations 
are  shown  schematically  in  Fig.  2.1,  and  in  the  as-built  construction 
drawings  of  Fig.  2.2.  The  pressure  data  were  to  provide  the  desired 
local  effects  and,  when  averaged,  the  net  force  on  the  cylinder.  The 
geometry  of  the  cylinder  was  suoh  that  the  gages  at  mid-length  would 
provide  the  loading  on  an  essentially  two-dimensional  body,  at  least 
during  the  initial  period  of  loadiiig  (see  Section  3.1). 
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Strain  gages  were  mounted  on  the  front  and  back  of  each  sensor 
support,  and  connected  in  series  electrically  to  cancel  out  bending 
strains.  Corresponding  sensor  bars  on  opposite  sides  of  the  cylinder 
were  averaged,  so  that  a total  of  two  instrumentation  channels  per 
cylinder  were  employed  for  strain  measurements.  The  strain  data  pro- 
vides a measure  of  the  net  force  acting  on  the  cylinder.  However,  it 
should  be  recognized  that  the  strain  data  as  such  does  not  necessarily 
give  a direct  measure  of  the  applied  blast  loading,  but  rather  the 
response  of  the  cylinder-sensor  system  to  this  loading.  In  effect 
the  cylinder  is  being  utilized  as  a dynamometer,  and  a suitable  analy- 
sis of  the  strain  data  is  required  in  order  to  determine  the  desired 
information. 

Soil  stabilization  was  provided  for  a distance  of  I50  ft  in  front 
of  the  cylinders  at  the  intended  U900-ft  location.  The  purpose  of  the 
stabilization  was  to  increase  the  probability  of  obtaining  a clean  pres- 
sure wave  at  that  location.  No  stabilization  was  provided  for  the 
cylinders  at  the  farther  location. 

Six  air-pressure  gages  were  placed  on  the  roof  and  sides  of  the 
UPSHOT -KNOTHOLE  Project  3-26.1  (Tests  of  the  Effects  on  POL  Installations, 
WT-756)  vertical  storage  tank.  The  location  of  these  gages  is  shown 
schematically  in  Fig.  2.6. 


2.2.2  Air-Pressure  Measurements 


All  air-pressure  - versus  - time  measurements  were  obtained  by 
the  use  of  Wiancko-type  gages;  a differential  inductance  bridge  actuated 
by  a pressure-sensitive  Bourdon  tube.  The  gage  output  was  fed  into 
modified  Webster-Chlcago  magnetic  tape  recorders.  The  circuitry  is  des- 
cribed as  a phase-modulated  system. 

The  pressure  gages  were  calibrated  statically  in  conjunction 
with  the  recording  system  just  prior  to  each  shot.  A regulated  air- 
pressure  system  was  used  for  positive  pressures  and  a vacuum  pump  for 
negative  pressures.  Accurately  known  steps  of  pressure  were  applied 
in  increments  of  10  per  cent  of  full-scale  deflection  for  each  gage. 

The  resulting  record  was  then  played  back  and  a calibration  established. 

The  accuracy  of  the  pressure  values  is  estimated  by  the  Ballistic 
Research  Laboratories  to  be  5 per  cent^  of  full-scale  readings.  The 
time  resolution  is  in  all  oases  within  one  millisecond.  Complete  details 
of  the  pressure  gage  installations  is  contained  in  the  final  report  on 
UPSHOT -KNOTHOLE  Project  5-28.1  (Structures  Instrumentation,  WT-758). 


2.2.5  Strain  Measurements 


Standard  SRlf  strain  gages  were  used  in  a four  active  am  bridge 
configuration  to  measure  axial  strain.  The  output  of  each  bridge  was 
fed  into  a Webster-Chicago  recording  system  through  a coupling  unit  and 


2/  This  is  not  felt  to  be  a realistic  error  estimate  insofar  as  the 
uncertainty  associated  with  information  obtained  from  the  pressure 
data  is  concerned.  In  fact,  ARF  believes  that  It  is  impossible  to 
assign  a rational  bound  to  the  probable  error  in  pressure  magnitudes. 
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recorded  on  magnetic  tape.  The  electrical  calibration  of  the  strain 
gages  was  accomplished  by  shunting  the  proper  aim  of  each  gage  installa- 
tion with  an  accurately  known  resistance  to  simulate  actual  strain. 

Complete  details  of  the  strain  gage  installations  are  contained 
in  the  report  on  UPSHOT -KNOTHOLE  Project  5*28.1  (Structures  Instrumenta- 
tion, WT-758). 

2.2.4  Instrument  Records  m 

I 

The  Ballistic  Research  Laboratories  handled  all  of  the  instru- 
mentation on  this  program.  The  output  of  the  pressure  and  strain  gages 
was  recorded  initially  on  magnetic  tape  and  later  played  back  onto 
oscillographic  paper.  The  resulting  records  had  nonlinear  ordinate 
scales  and  were  undesirable  for  purposes  of  interpretation  and  compari- 
son. Thus,  all  the  records  were  converted  into  linear  form. 

The  BRL  reduced,  calibrated,  and  plotted  to  linear  scales  all 
of  the  instrument  records.  The  records,  as  presented  to  ARF,  were  in 
the  form  of  plots  made  up  of  the  points  at  which  the  records  were  read. 

The  ARP  was  responsible  for  fairing  curves  through  these  points.  Hie 
BRL  also  submitted  tabulated  listings  of  the  points,  as  well  as  copies 
of  the  original  playbacks. 

The  numerical  data-reduction  scheme  employed  on  the  strain 
records  required  that  the  data  be  given  at  equal  time  intervals. 

Accordingly,  the  BRL  made  a separate  linearization  of  selected  strain 
records  at  equal  time  intervals  of  about  k ms,  and  supplied  these 
results  on  IBM  punched  cards. 

2.5  LOCATION  OF  TEST  STRUCTURES 

The  location  of  the  cylinders  at  the  test  site  is  shown  in 
Fig.  2.7.  Ground  range  from  intended  and  actual  ground  zero  is  given 
in  Chapter  4,  Table  4.1.  The  measured  overpressure  levels  are  eiIso 
listed  in  that  table.  The  orientation  effects  resulting  from  the  bomb- 
ing error  are  considered  to  be  negligible.  The  test,  as  conducted, 
did  not  significantly  deviate  from  the  test  specifications  given  in 
the  final  pretest  report  (Air  Force  Planning  Program,  Final  Pretest 
Report,  Part  VI,  Tests  on  Horizontal  Cylinders,  with  the  exception 
that  pretest  planning  anticipated  participation  in  only  one  shot. 
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a.  Typical  Cylinder  Shoving  Pressure  and  Strain  Agje  Locations 


Fig.  2.2  As-Bullt  Construction  Drawings 
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Schedule  for 


2.3  1-1/2  ft 
Diameter  Cylinder 
3.3b 


Fig.  2.4  End  View  of  Typical 
6 ft  Diameter  Cylinder 
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PRETEST  CaHSIIERATICMS 


3.1  SETEHTTriW  (g  TEST  HEMS 

The  xxamber  aaai  location  of  the  cylinders  were  chosen  so  as  to 
achieve  the  test  objectives  with  the  possible  number  of  cylin- 
ders. In  the  test,  the  following  items  were  taken  into 

consideration. 

Shock  tube  eocperiments  at  Mnceton  University  (Bleakney,  W. 
Diffraction  of  Shock  Waves  Around  Obstacles.  Technical  Report  II-3) 
IxkLlcate  that  the  character  of  the  loading  on  the  top  of  two-dimenr* 
sional  cylinders  of  square  cross  section  is  dependent  on  the  strength 
of  the  shock.  This  is  apparently  due  to  the  dependence  of  the  formation 
of  vortices  with  shock  strength,  and  it  is  quite  possible  that  a similar 
effect  occurs  in  the  case  of  circular  cylinders.  This  effect  was 
investigated  in  the  test  by  jdacing  two  cylinders  of  the  same  size  in 
overpressure  regions  differing  by  about  a factor  of  2 (Cylinders  3 •3a 
aai  3.3c).  It  was  felt  that  any  vorteac  effects  dependent  on  shock 
strength  should  be  discernible  for  this  variation  in  pressxn'e. 

Shock  tube  eacperimeints  indicate  that  the  time  details  of  the 
inarHwg  dtalng  the  diffraction  phase  scale  linearly  with  the  absolute 
size  of  the  geoaetry  under  consideration.  This  dependence  of  size  was 
investigated  by  placing  two  cylinders  of  the  same  diemeter-to-length- 
to»grouai  proocimity  ratio,  but  differing  in  absolute  size  by  a factor 
of  4,  in  the  same  overpressure  region  (Cylinders  3.3a  and  3.3b).  A 
scale  ratio  of  4-to-l  was  felt  to  bo  sufficient  in  order  to  point  out 
this  effect  of  size. 

In  addition  to  the  scaling  of  time  details  during  the  diffraction 
period,  the  loading  during  the  drag  phase  is  also  dependent  upon  the 
absolute  size  of  the  cylinder.  The  total  horizontal  steady-state  drag 
coefficient  for  cylinders  is  known  to  vary  between  1 and  1/3  in  the 
range  of  Reynolds  numbers  (based  on  diameter)  between  ICP  and  7 x 10^, 
then  remain  essentially  constant  at  1/3  for  larger  Reynolds  numbers 
(see  AerQdvTtf7ir<^r.  by  S.  F.  Hoemer).  Under  transient  flow  condi- 
tions such  as  in  the  field,  the  Reynolds  number  itself  will  vary  with 
time.  For  the  test  cylinders,  the  Refolds  nuiober  is  supercritical 

(i.e.,  greater  than  7 x 105)  and  will,  of  course,  vanish  vdien 
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the  porticle  velocity  goes  to  aero.  There  is  considerable  interest  in 
the  behavior  of  the  drag  loading  under  these  circumstances.  However, 
there  is  some  question  as  to  whether  the  pressure  gage  instrumentation 
is  adequate  to  detenaine  this  effect  of  siae.  That  is,  under  test 
conditions  the  drag  pressures  are  expected  to  be  smll  - of  the  order 
of  uncertainty  in  the  pressure  datum  itself.  Therefore,  infome-tion 
of  this  nature,  \iiile  highly  desirable  is  considered  to  be  of  bonus 
value. 

Tests  conducted  in  the  i'VHF  4-in.  shock  tube  showed  that  the 
pressure  on  the  xaaderside  of  a cylinder  depends  markedly  on  the 
OToxindty  of  the  cylinder  to  a reflecting  surface  such  as  the  ground 
(Air  Force  Planning  Program,  Final  Pretest  Pteport,  Part  I,  Oil  Storage 
and  EouiTTTK^nt.  Teats') . Jt  is  believed  that  a systematic  study  is  in 
order  to  determine  quantitatively  both  the  magnitude  of  this  pressure 
build  up  and  the  scaling  relationships  which  hold  between  full-scale 
and  model  testa.  Three  cylinders  of  the  same  size  but  at  different 
distances  above  the  ground,  subject  to  the  same  loading,  \wre  to  pro- 
vide the  desired  information  (Cylinders  3.3c,  3.3d,  and  3.3e). 

It  is  felt  that  a reasonable  comparison  of  two-  and  three- 
dimensional  loadings  could  be  had  by  cor^nring  the  net  loads  as 
deterzained  by  the  pressure  data  (two-dimensional  loading)  and  ly  the 
strain  data  (three-dimensional  loading).  The  ratio  of  cylinder  length 
to  diameter  is  believed  to  be  sufficiently  great  (L/D  = 3.33)  so  that 
the  pressure  gages  located  at  mid-length  will  read  essentially  two- 
dimensional  data,  at  least  for  the  initial  period  of  the  loading.  The 
pressure  gages  should  represent  two-dimensional  loading  exactly  for  the 
first  3r/U  units  of  time  (r  s radius  of  cylinder,  U ■ velocity  of  the 
shock  front) , since  it  woiild  take  a disturbance  signal  at  least  this 
length  of  time  to  travel  from  one  end  of  the  cylinder  to  the  center, 
llbreover,  it  is  believed  that  such  signals  are  sufficiently  weak,  so 
that  the  observations  made  at  mid-length  reptresent  the  two-dimensional 
load  on  the  cylinder  for  the  entire  diiration  of  the  loading.  The  data 
obtained  from  the  strain  gages,  on  the  other  hand,  are  indicative  of 
the  net  loading  on  a finite  cylinder  and,  as  such,  represent  the  thre»- 
dimensional  loading  case.  All  of  the  test  cylinders  may  be  said  to 
contribute  to  this  objective. 

3.2  MAI) 

Theories  concerning  the  influence  of  ground  prcodmity  on  the  blast 
loading  on  cylindrical  shapes  are  practically  nonexistent.  For  a 
circular  cylinder  in  free  space  (i.e.,  infinitely  removed  from  the 
ground)  the  situation,  is  more  satisfactory.  Prediction  methods  for  this 
case  are  discussed  in  detail  in  the  final  pretest  report  (Air  Force 
Planning  Program,  Final  Pretest  Report,  Part  VI,  Tests  on  Horizontal 
CvUnders^ , and  the  results  are  presented  in  Chapter  5 of  this  report 
(see  Fig.  5.13) . The  tall  chimney  included  in  Operation  JAHGLE  (Air 
Force  Structures  Program  3.3»  Final  Report,  Vol.  l)  affords  a check 
between  measured  and  predicted  loads  for  this  case. 
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The  pretest  work  on  load  predictions  was  concerned  priwarily  with 
an  estimate  of  the  peak  pressure  on  the  underside  of  the  front  portion 
of  a cylinder  as  a function  of  ground  proximity.  No  attempt  was  made 
to  predict  the  subsequent  pressure  distribution  as  the  shock  spills 
over  the  rear  of  the  cylinder.  This  work  was  based  primarily  on  the 
results  of  shock  tube  studies  (see  Bleakney,  W.  The  Diffraction  of 
Shock  Waves  Around  Obstacles,  Technical  Report  II-5,  Princeton  Univer- 
sity and  Air  Force  Planning  Program,  Final  Pretest  deport.  Part  I, 

Oil  Storage  and  Equipment  Tests) . The  limiting  case  for  which  the 
cylinder  touches  the  ground  has  been  solved  analytically  from  an  approxi- 
mate point  of  view.  This  condition  was  assumed  to  be  equivalent  to  the 
pressure  twlld  up  occxirring  \^en  a shock  impinges  upon  a re-entrcnt 
corner  (Air  Force  Planning  Program,  Final  Pretest  Report,  Part  V, 

Regular  Reflection  on  Cubical  Structures).  The  other  limiting  case  is, 
of  course,  for  the  cylinder  in  free  space.  Inasmuch  as  these  predic- 
tions are  rather  crude  and  no  comparison  with  test  results  proved  possi- 
ble, the  work  is  not  repeated  here. 

5.3  NET  FORCE  MEASUREMENTS 

3.3»1  Sensor  Design 

Effort  was  expended  during  the  planning  phase  of  the  program 
to  devise  a method  for  direct  measurement  of  the  net  applied  blast 
forces  acting  on  the  test  cylinders.  Ideally,  a direct  measurement 
system  would  be  one  whose  output  is  either  proportional  to  the  net 
blast  loads  acting  over  the  entire  structure,  or  whose  output  can  be 
converted  to  net  loads  in  a straightforward  fashion.  This  requirement 
for  a so-called  net  force  measurement  system  is  meant  to  exclude  pres- 
sure gage  instrumentation,  at  least  in  the  present  case.  While  there 
are  certain  types  of  structures  for  which  net  forces  (over  the  entire 
structure)  can  be  adequately  determined  by  means  of  combining  and 
averaging  individual  pressure  gage  data,  there  are  many  structures  of 
interest  where  such  averaging  is  either  subject  to  considerable  uncer- 
tainty or  is  totally  impractical.  There  are  also  cases  where  the 
presence  of  pressure  gages  may  sufficiently  alter  the  flow  conditions 
around  the  test  item,  so  as  to  rule  out  this  means  of  force  measurement. 
In  the  present  case,  where  an  averaging  procedure  was  possible,  it  was 
not  practical  to  provide  the  necessary  additional  pressure  gage  instru- 
mentation. Thus,  an  alternative  force  measurement  system  had  to  be 
devised. 

The  force  measurement  or  sensor  system  finally  adopted  con- 
sisted of  supporting  the  cylinders  on  either  end  as  shown  in  Fig.  2.2, 
and  measuring  the  forces  in  these  supports  by  means  of  resistance-type 
strain  gages.  While  other  supporting  schemes  were  considered  (e.g., 
a bsilllstlc  pendulum  arrangement  utilizing  accelerometer  output)  the 
present  scheme  was  resorted  to,  since  an  identical  system  was  utilized 
in  connection  with  the  UPSHOT -KNOTHOLE  Project  3.U  (Tests  on  the  Loading 
of  Truss  Systems  Common  to  Open-Framed  Structures,  WT-725 ) . In  this 
latter  program  the  end- support  scheme  appeared  to  be  the  only  practical 
approach . 
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In  the  present  setup,  the  test  cylinder  was  utilized  in  the 
sense  of  a dynamometer.  That  is,  the  cylinder  and  sensor  arrangement 
as  a whole  constituted  the  gaging  device.  The  output  of  the  strain 
gage  system  (i.e.,  the  force  in  the  sensor  bars)  was  thus  dependent 
upon  the  dynamic  characteristics  of  the  cylinder-sensor  system  as  well 
as  upon  the  applied  blast  force.  In  order  for  the  strain  data  to  be 
utilized  as  a direct  measure  of  the  applied  blast  loadings,  it  is  neces- 
sary that  the  frequency  response  of  the  cylinder  be  considerably  greater 
than  that  of  the  sensor  supports  alone  and,  in  turn,  that  the  funda- 
mental period  of  the  system  as  a whole  be  sufficiently  less  than  the 
significant  time-details  of  the  applied  blast  loading.  In  this  case 
the  response  is  essentially  that  of  a single-degree-of -freedom  system 
(i.e.,  a rigid  cylinder  on  flexible  supports)  and,  with  suitable  damp- 
ing, the  strain  data  could  be  interpreted  much  as  axe  the  pressure  data  - 
simply  by  inspection.  Unfortunately,  it  did  not  appear  possible  to 
satisfy  these  conditions.  While  the  rigidity  of  the  cylinders  could 
have  been  increased  considerably,  the  frequency  response  of  the  sensor 
bars  was  limited  in  two  ways.  First,  was  the  requirement  that  a suf- 
ficiently great  strain  be  produced  to  yield  a reliably  measurable 
signal  from  the  strain  gages  throughout  most  of  the  loading  period. 
Second,  the  sensor  arrangement  was  designed  on  the  basis  of  there  being 
no  provision  for  pretest  field  '.alibration  of  the  entire  setup  (and 
there  was  none).  This  necessitated  a rather  longer  sensor  bar  than 
might  otherwise  have  been  utilized  (and,  hence,  a lower  frequency), 
since  it  was  Imperative  that  the  measured  surface  strain  be  an  inade- 
quate measure  of  the  actual  force  in  the  bar.^ 

Details  of  the  sensor  bar  construction  are  shown  in  the  as- 
built  drawings  of  Figs.  2.2  and  in  the  photograph  of  Fig.  2.4.  Briefly, 
the  end  supports  consisted  of  a rigid  right -angle  frame  of  welded  con- 
struction whose  vertex  was  attached  at  the  coclal  center  of  the  cylinder. 
On  the  6-ft  diameter  cylinders,  the  bars  were  2 in.  square  and  approxi- 
mately 22  in.  long  between  welds  (i.e.,  the  unsupported  sensor  bar 
length);  on  the  quarter-scale  cylinder  (3-3h),  the  bars  were  1 in. 
square  and  approximately  11  in.  long.  High-strength  steel  was  used  for 
all  sensor  bars.  The  geometry  of  the  supports  was  such  that  each  sensor 
bar  could  be  considered  as  a two-force  member.^  This  would,  of  course, 
be  the  case  if  the  frame  were  Ideally  pin-connected,  but  a rigid  frame 
was  more  desirable  from  the  construction  viewpoint. 

"y  That  is,  the  sensor  bar  had  to  be  sufficiently  long  so  that  the 

stress  distribution  at  the  cross  section  being  gaged  was  essentially 
uniform. 

2/  An  analysis  of  the  maximum  error  Introduced  by  the  assumption  of 
pin-ends  was  carried  out  in  the  final  pretest  report  (Air  Force 
Planning  Program,  Final  Pretest  Report,  Part  VI,  Tests  on  Horizontal 
Cylinders) . The  error  in  the  magnitude  of  the  force  was  found  to 
be  less  than  6 per  cent,  and  about  1.5  degrees  in  direction  relative 
to  the  centerline  of  the  bars. 
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5.3-2  Interpretation  of  Strain  MeaBurements 


The  final  sensor  design  was  such  that  the  measured  forces  In 
the  cylinder  support  would  most  likely  be  Influenced  by  the  dynamic 
characteristics  of  the  entire  test  structure  and,  thus  would  net  be 
a direct  measure  of  the  applied  blast  force.  This  situation  necessi- 
tated a data-reduction  system  that  was  capable  of  taking  this  influence 
into  account  and  providing  the  desired  net  applied  forces.  In  essence, 
the  problem  here  is  the  inverse  of  the  conventional  dynamic  response 
problem.  That  is,  the  usual  problem  is  "given  the  dynamic  system  and 
the  input,  find  the  output"  — here  the  problem  is  "given  the  dynamic 
system  and  the  output,  find  the  input." 

The  solution  to  the  inverse  problem  is  generally  the  more 
ccanplicated  but  in  principle  can  be  achieved,  at  least  for  linear 
dynamic  systems,  provided  the  response  of  the  system  to  a given  input 
is  known.  However,  the  present  problem  is  made  more  difficult  by  the 
fact  that  the  characteristics  of  the  dynamic  system  itself  are  not 
known  with  any  certainty.  That  is,  the  system  is  sufficiently  compli- 
cated that  the  response  (l.e.,  the  force  in  the  sensor  bars)  to  a known 
input  cannot  be  computed  with  sufficient  accuracy  to  Justify  a solution 
to  the  inverse  problem  by  standard  methods. 

A data-reduction  scheme  referred  to  as "transient  analysis"  was 
developed  during  the  course  of  the  pretest  work  and  was  believed  to  be 
adequate  for  the  needs  of  this  program.  In  principle,  the  method  of 
transient  analysis  is  applicable  to  the  determination  of  the  time  vari- 
ation of  an  input  function  (the  magnitude  of  the  input  can  be  deter- 
mined only  within  an  arbitrary  scale  factor)  from  a known  output  function, 
even  when  the  equations  of  motion  of  the  system  are  not  known  in  detail. 
The  method  is  restricted  to  dynamic  systems  characterized  by  linear 
second-order  differential  equations  with  constant  coefficients,  and 
input  functions  which  can  adequately  be  represented  by  finite  sums  of 
exponential  functions  (i.e.,  "e"  functions),  all  of  whose  coefficients 
and  exponents  are  real  numbers.  In  effect,  this  last  requirement  limits 
the  input  functions  to  those  having  no  harmonic  components,  which  seemed 
a reasonable  enough  assumption  in  the  case  of  blast  loading  functions. 

A further  restriction  is  that  the  forcing  functions  associated  with 
the  various  modes  of  the  system  be  either  zero  or  timewise  proportional. 

The  restrictive  assumptions  which  permit  the  valid  use  of  tran- 
sient analysis  as  a technique  for  data  reduction  in  the  present  applica- 
tion were  considered  to  be  satisfied  during  the  pretest  planning  period. 
However,  as  discussed  in  Chapter  3,  the  application  of  transient  analysis 
to  the  strain  data  was  unsuccessful.  The  analysis  is  developed  in 
Appendix  A of  this  report,  which  also  contains  a critical  discussion  of 
the  method. 
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Relationship  Between  Strains  and  Forces 
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The  relationship  between  strain  and  total  force  in  the  sensor 
bars  per  unit  projected  area  of  the  cylinder  is  given  by 

S=(A^A)E€  (5.1) 

5 C 

where, 


S 

A 

s 

A 

c 

E 

€ 


unit  force  in  one  sensor  bar,  psi 

2 

cross-sectional  area  of  the  bar,  in. 


projected  area  of  the  cylinder  along 
a diameter,  in.^ 

Young's  modulus  for  the  sensor  bars,  psi 
strain  measured  in  one  sensor  bar,  in. /in. 


Assuming  pin  connections,  ‘•.he  horizontal  and  vertical  components 
of  to'tal  load  per  unit  projected  area  acting  on  the  cylinders  is  given 
in  terms  of  the  force  in  the  sensor  bars  as : 


F 


h 


F 

V 


+ Sg) 


(5.2) 


The  unit  horizontal  force,  Fj^,  is  taken  to  be  positive  acting  down- 
stream, and  the  unit  vertical  force,  F , positive  acting  upward.  The 
subscripts  "1"  and  "2"  refer  to  the  forces  in  the  upstream  and  down- 
stream sensor  bars,  respectively,  which  are  assumed  positive  in  tension. 
In  terms  of  'the  measured  strains,  F^^  and  F^  are  given  numerically 
as. 


F 


h 


F 


V 


9.82  X 10“^(e^  - eg),  psi 

> 

9.82  X 10‘^(«^  + eg),  psi 


for  the  6 -ft  diameter 
cylinders,  5.5®>  c,  d, 
and  e 


(5.5a) 


- 59.5  X 10"5(€^  - 6^),  psi 

- 59.5  X 10"^ (e^  eg),  psi 

where  the  strains  are  given  in  ji  ln./in.(»  10"*^  in. /in.),  and  a value 
of  E = 50  X 10®  psi  has  been  used  for  Young's  modulus. 
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for  the  1.5-ft  diameter 
cylinder,  5.5b 


(5.5b) 
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CHAPTER  4 


EXPERIMENTAL  RESULTS 


k.l  FIELD  CONDITIONS 

The  specifications  for  the  flow  conditions,  as  established  during 
the  pretest  planning  period,  were  not  substantially  violated  during 
the  test.  The  cylinders  were  located  in  the  Mach  reflection  region 
in  both  Shots  9 and  10,  and  the  incident  shock  waves  were  clean  and 
well -developed  in  each  case.  The  error  in  ground  zero  on  Shot  9 
changed  the  range  and  orientations  of  the  cylinders,  and  resulted  in 
somewhat  different  overpressures  than  were  intended.  (Shot  10  was  not 
anticipated  during  the  planning  phase.)  The  field  conditions  are 
summarized  in  Table  4.1,  and  the  blast  wave  parameters  are  given  in 
Table  4.2. 

4.2  PRESSURE  GAGE  DATA 

Of  the  thirty  pressure  gages  used  in  Shot  9>  four  did  not  yield 
pressure-time  observations,  whereas  of  the  thirty  gages  used  in  Shot  10 
only  one  failed  to  yield  any  observations.  The  records  lost  on  Shot  9 
were  frcm  Gages  P2  and  P5  of  Cylinder  5- 5°  and  from  Gage  P5  on  Cylin- 
der 3.3d  and  3.3e.  The  only  record  lost  on  Shot  10  was  from  Gage  P5 
of  Cylinder  3.3e.  (See  Fig.  2.1  for  pressure  gage  locations.) 

In  general,  the  original  pressure-time  observations  which  were 
obtained  appeared  to  be  quite  good  on  first  inspection;  that  is,  the 
noise  level  was  reasonably  low  and  the  records  were  relatively  free 
from  extraneous  pips,  etc.  These  data  are  considered  in  more  detail 
in  Chapter  5.  Typical  pressure-time  records  are  shown  in  Fig. 

4.3  STRAIN  GAGE  DATA 

Of  the  twenty  strain  gages  (ten  each  on  Shots  9 aud  10),  all  but 
one  provided  records  of  some  sort  in  each  shot.  The  records  were  played 
back  by  BRL  at  two  speeds;  a slow  speed  which  in  most  cases  gave  infor- 
mation during  and  beyond  the  positive  leading  period,  and  a fast  speed 
which  indicated  greater  detail  of  the  response  during  the  first  200  ms. 
or  so  after  shock  arrival.  The  extent  of  these  two  playbacks  is 
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TABLE  4.1  - Summary  of  Field  Conditions 


Ground 

Orientation 

Actual 

Test 

Shot 

Range 

(Degree  frcmi 

Overpressure 

Duration* 

Item 

No. 

(ft) 

normal) 

(psig) 

(sec. ) 

Intended 

Actual 

Intended 

Actual 

Intended 

Actual 

9 

4900 

4780 

0 

9-1/2 

8 

6.7 

0.90 

3.3a 

4900 

4750 

0 

1 

- 

3.5 

0.90 

9 

4900 

4780 

9-1/2 

8 

6.7 

0.90 

3.3b 

4900 

4750 

1 

- 

3.5 

0.90 

9 

6375 

6280 

7-1/2 

5.2 

4.6 

0.95 

3.3c 

6375 

6290 

0 

1 

- 

2.0 

0.95 

9 

6375 

6280 

0 

7-1/2 

5.2 

4.6 

0.95 

3.3«i 

6375 

6290 

0 

1 

2.0 

0.95 

9 

6375 

6280 

0 

7-1/2 

5.2 

4.6 

0.95 

3.3e 

6375 

6290 

0 

1 

- 

2.0 

0.95 

Ambient  preshot  conditions:  Shot  9 air  pressure  = IJ.l  psi; 

air  temperature  = l6.7“  C. 

Shot  10  air  pressure  = 15.2  pslj 
air  temperature  = l4.8°  C. 

* 

Duration  of  the  Positive  Phase  of  the  Blast  Wave. 

Data  obtained  frcai  Summary  Report  of  the  Technical  Director,  WT-782. 


34 

HeSTKteTEP  BAIA 

UNCLASSIFIED 

s.  ^ 


table  4*2  - Sussaary  of  Blast  Wave  Farazoetars 


Shot  9 

Shot  10 

Parameter 

Symbol 

Units 

5.5^ 

3.5c,d 
and  e 

3.5c,d 
and  e 

Ambient  Pressure 

psi 

13.2 

13.2 

13.2 

13.2 

Ambient  Temperature 

’’o 

•P 

62 

62 

59 

59 

Ambient  Sound  Velocity 

''o 

fps 

1121 

1121 

UI8 

lllB 

Ambient  Density 

Po 

Ib/ft^ 

0.0683 

0.0685 

0.0685 

0.0683 

Peak  Overpressure 

pel 

6.7 

4.6 

3.5 

2.0 

Shock  Strength 

i 

- 

1.51 

J-.35 

1.26 

1.15 

Reflected  Pressure 

Pr 

psl 

16.1 

10.6 

7.67 

4.09 

Reflected  Pressure 

Ratio 

Pi./Po(0) 

- 

2.40 

2.30 

2.23 

2.06 

Duration  of  Positive 
Phase 

^0 

sec 

0.90 

0.95 

0.90 

0.95 

Density  Behind  Shock 

Pi 

Ib/ft^ 

0.0916 

0.0841 

0.0800 

0.0752 

Temi)erature  Behind 

Shock 

•f 

127 

107 

95 

78 

Particle  Velocity 

Behind  Shock 

u 

fps 

542 

24? 

188 

112 

Shock  Front  Velocity 

u 

fps 

1341 

1280 

1254 

1185 

Peak  Drag  Pressure 

Pd(o) 

psi 

1.13 

0.53 

0.30 

0.11 

Viscosity 

fx  X 10^ 

Ib/ft  sec. 

133 

129 

127 

125 

Reynolds  Number 
(based  on  Diameter) 

- 

14^2^ 

^'5^2 

9.65 

T^Og^ 

4.02 

Notes:  Ambient  conditions,  overpressure  and  duration  obtained  from  Sunmary  Report 
of  the  Technical  Director,  Wr-782. 

Other  quantities  ccmputed  from  the  conventional  Ranklne-Hugonlot  relationships 
applicable  across  the  shock  front. 
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TABLE  4 >3  - Smnnary  of  Strain  Gage  Data 


Strain 

1 Remarks  (see  notes  at 

end  of  table) 

Record 

Shot  9 

Shot  10 

3.3*S1 

s > t^;  f (NL)  fif  180  ms 
i ECg 

s > t^;  f a 1U5  me 

. BCg 

3.5«82 

b(l)  > t^;  f 2;  lUO  ms 
cr  ECg 

S < t^;  f 2f 300  ms 

“l-*=2 

3.3b61 

No  record 

B < t^;  f cr  190  ms 

3.3bS2 

8 > t : t — 200  ms 
o' 

ECg 

*(L)  > t^:  f 2f  210  ms 

f and  8 do  not  agree 
Initially 

ECj^Sf  BCg 

e > t^;  f(NL)  Cr  350  ms 

BCi^EC2 

5.5cS2 

s < t^;  f 520  ms 

ECj^ar  ECg 

No  record 

5.3dSl 

B srt^;  f 2r260  ms 

BC^  4 ECg 

8 > t^i  f 260  ms 

®Cl=*=2 

5.3dS2 

s < t^;  f Cr  270  ms 

EC  cr  EC_ 

1 2 

* > t ; f 260  ms 

0 

3.5eSl 

s > t^;  f of  270  ms 

ECi  4 Kg 

s < t^;  f 2r380  ms 

EC^-  ECg 

3.5eS2 

a < t^;  f Cf 300  ms 

ECj^  - ECg 

8 > t ; f 2f 2U0  ms 

0 

ECj^  . BCg 

Note*:  1.  See  Fig.  2.1  for  explanation  of  strain  gage  code. 

2.  The  ten)  "*"  refer*  to  elov  record  playback.  "*(L)"  Indicates  a linear- 
ized plot,  otherwise  ”*"  plots  not  linearized,  "s  > t^"  Indicates  "s"  plot 
extends  beyond  duration  of  positive  loading  period,  otherwise  "s  < t ". 

3.  "f"  refer*  to  fast  record  playback.  "f(NL)"  Indicate*  "f"  plot  not  ° 
linearized,  otherwise  all  "f"  plots  linearized,  "f  flf  l80  ms"  Indicates 
"f"  plot  extends  to  about  l8o  ms  after  shock  arrival. 

I*.  BC.,  BCg  refer  to  pretest  electrical  calibrations  (see  text). 


indicated  in  Table  idilch  BunmaplzeB  all  the  Btraln  gage  data. 

Only  two  of  the  alow  playbackB  were  linearized  (5.5aS2,  Shot  9> 

5.3bS2,  Shot  10),  whereaB  all  but  two  of  the  faet  playbackB  were 
linearized. 

The  entrieB  in  Table  U.3  labeled  EC^  and  EC2  refer  to  electrical 
calibration  marka;  ECi  1b  a field  calibration  Juat  prior  to  ahot  time, 
and  ECg  is  an  earlier  calibration  taken  Bome  tine  before  the  teata. 

No  poBtahot  callbratlona  were  made.  Ab  1b  Been,  theee  marke  are 
generally  only  approximately  equal,  and  in  thoae  caaea  marked  "not 
equal,"  the  difference  waa  of  the  order  of  the  atrain  calibration 
Jumps.  The  llnearlzationa  were  made  on  BRL'a  beat  eatimate  of  the 
proper  calibration;  where  BCi  differed  appreciably  from  BC2,  the  poa- 
aible  error  waa  not  eatlmated  by  BRL.  Aa  far  aa  could  be  told,  the 
baaellnea  appeared  to  be  good  in  that  no  drift  from  the  preahot  poal- 
tion  could  be  noted  by  vlaual  Inapectlon. 

Aa  Indicated  in  Table  U.J,  the  magnitude  of  Initial  peaka  on 
record  5.3bS2,  Shot  10,  do  not  agree  between  the  faat  and  alow  play- 
backa.  This  was  the  only  case  in  -vrtiich  a discrepancy  between  a linear- 
ized plot  and  the  original  playback  could  not  be  resolved.  However, 
not  evident  in  the  table  is  the  fact  that  many  similar  diacrepancies 
arose.  For  example,  three  linearized  plots  were  originally  discarded, 
since  they  appeared  to  provide  data  for  less  than  35  mOj  until  it  was 
noted  that  the  time  scales  were  in  error  by  a factor  of  10.  Additional 
discussion  of  the  reliability  of  the  strain  da+a  is  given  in  Chapter  5- 
Typical  strain  records  are  shown  in  Fig.  5«28. 
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CHAPTER  5 


DISCUSSION  OF  RESULTS 


5.1  PRESSURE  MEASUREMEafTS 

5.1.1  Reliability  of  Data 

Spot-checking  the  original  record  playbacks  with  the  BRL 
linearizations  showed  some  discrepancies.  These  consisted  mainly  of 
extraneous  pips  in  the  linearizations  and  were  easily  corrected  by 
comparing  them  with  the  original  playbacks.  However,  several  not  so 
obvious  irregularities  were  noted,  aud  it  is  probable  that  some  of 
these  were  carried  through  the  analysis.  A typical  example  of  this 
type  of  discrepancy  is  seen  by  comparing  Figs. 5.2a  and  b.  Figure  5.2a 
is  a partial  tracing  of  the  original  playback  on  Record  5.5eP5,  Shot  10; 
Fig.  5.2b  is  the  linearized  plot  of  the  same  record.  As  can  be  seen, 
the  magnitudes  of  the  first  and  second  maximums  and  the  first  minimum 
are  not  consistent  between  these  records;  after  about  15  ms,  however, 
the  two  compare  favorably. 

A more  serious  criticism  of  the  records  is  that  the  magnitudes 
of  the  pressures,  in  general,  are  lower  than  expected.  The  difference 
is  apparent  in  the  few  instances  where  lower  bounds  can  be  established 
for  the  pressures.  For  example,  at  the  position  of  Gage  P5  (the  head- 
on  gage  position)  the  pressure  should  rise  instantaneously  to  the  normal 
reflected  pressure,  p^.,  and  then  decay  to  its  pseudo- steady-state 
value.  As  seen  in  Fig.  5-5^  this  gage  on  Cylinders  5.3*  and  3. 3c, 

Shot  10,  recorded  peak  pressures  of  between  1.6pg.(0)  and  1.7p^(0),  where- 
as the  expected  values  (listed  in  Table  4.2)  were  about  2.2'p^{0)  and 
2.1pg^0),  respectively.  On  Cylinder  3*3®,  Shot  9,  Gage  P3  gave  a value 
of  Pr/P{T(0)  » 2.1,  which  is  12  per  cent  lower  than  the  expected  value 
of  2.6. 

At  the  head-on  gage  position,  the  pseudo-steady-state  value  is 
always  greater  than  the  free-stream  pressure,  Pa(t),  since  the  drag 
coefficient  is  positive  at  that  point.  The  free-stream  pressure-time 
curves  (obtained  from  Summary  Report  of  the  Technical  Director,  WT-T82) 
are  presented  in  Fig.  5-9  for  about  the  first  12 5 ms  after  shock 
arrival.  For  convenience,  these  curves  have  been  normalized  in 
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terms  of  peak  side-on  pressure  and  r/U  time  units,  where  r = 5 f t 
is  the  radius  of  the  larger  cylinders.  A comparison  of  Figs.  5*5  and 
5.9  shows  that  the  pressures  recorded  by  Gage  P5  on  Cylinders  5- 3a, 

Shot  9,  and  5 *30,  Shot  10,  are  slightly  greater  than  free-stream  pres- 
sure at  the  time  t • 50  r/U  as  expected,  while  the  same  gage  on 
Cylinder  5 •5a,  Shot  10,  reads  just  slightly  less  than  free-streeun 
pressure. 

The  initial  pressure  at  the  positions  of  Gages  P2  and  P4  (see 
Fig.  2.1)  should  also  be  of  the  order  of  the  reflected  pressure,  p_, 
whereas  the  pseudo-steady  state  value  of  the  pressure  at  Gage  P2  should 
be  approximately  p^(t),  since  the  drag  coefficient  at  that  position 
is  approximately  zero.  As  seen  in  Fig.  Gage  P2  indicates  a peak 

value  less  than  1.7p_(0)  for  Cylinder  5*3a,  Shots  9 and  10,  and  for 
Cylinder  5-5c,  Shot  &.  At  a time  of  50  r/U  units.  Gage  P2  reads  less 
than  free-stream  pressure  in  both  shots;  the  same  gage  on  Cylinder  5«5c, 
Shot  10,  is  approximately  equal  to  free-stream  pressure  at  this  time. 

Gage  P4  on  Cylinder  5-5a,  Shot  9 (Fig.  5*6),  appears  reasonable  in 
Shot  9 and  probably  somewhat  low  in  Shot  10.  However,  this  gage  on 
Cylinder  5 *30  reads  a peak  value  less  than  the  incident  free-stream 
pressure  in  both  shots.  While  it  is  possible  that  the  gage  system  did 
not  respond  to  the  initial  peak  pressures,  this  seems  unlikely  and 
would  not  account  for  the  lower  pseudo- steady- state  values  noted  above. 

The  recording  systems  were  electrically  calibrated  twice: 
once  some  time  before  the  test  (mark  EC2  in  Fig.  5»2a)  and  again  just 
prior  to  the  test  (mark  EC.  in  Fig.  5*2a).  These  calibration  steps 
did  not  match  in  approximately  half  of  the  records;  however,  the 
majority  of  the  deviations  were  small.  Thus,  the  same  uncertainties 
found  in  the  strain  data  calibrations  noted  in  Chapter  k are  present 
here. 

The  original  data  appeared  to  contain  also  some  shifts  in  base- 
line. The  combination  of  the  baseline  shifts  and  calibration  errors 
affect  the  magnitudes  of  the  pressure  records,  and  since  the  scales 
were  nonlinear,  these  errors  can  affect,  to  a small  degree,  the  shape 
of  the  records.  The  time  scales,  however,  seem  to  be  quite  reliable. 

There  is  some  question  of  data  reliability  regarding  overpres- 
sures, since  the  actual  overpressures  were  lower  than  the  intended 
ovei^jressures  (see  Table  U.l).  This  resulted  in  pressure  deflections  on 
the  recorder  which  were  less  than  full  scale,  and  hence,  the  recorder 
and  calibration  errors  were  effectively  enlarged.  BRL  estimated  the 
accuracy  of  the  pressure  data  to  be  5 per  cent  of  full-scale  deflection.^ 

To  facilitate  analysis  of  the  relatively  large  number  of  pressure- 
time  records,  the  linearized  form  of  the  data  was  replaced  by  a "simpli- 
fied" pressure-time  curve.  Typical  examples  are  shown  in  Fig.  5«1* 

A series  of  straight  lines  were  drawn  through  the  linearized  data, 
conserving  the  essential  features  of  the  curves  and  neglecting  the 
oscillations  which  were  due  to  the  dynamic  effects  of  the  recording 
system.  The  accuracy  of  this  procedure,  of  course,  depends  upon  the 
analyst's  understanding  of  the  basic  flow  phenomenon  as  well  as  the 


^ See  footnote  on  page  19, 
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dynamic  recording  system.  These  simplified  curves  vere  plotted  in 
a normalized  form,  using  the  variable  p/Pjj.(0)  for  the  unit  of  pres- 
sure and  r/U  for  the  unit  of  time,  where  p is  measured  pressure, 
p (O)  is  the  peak  incident  overpressure,  r is  the  raiius  of  the 
c^inder,  and  U is  the  velocity  of  the  shock  wave.  This  final  form 
of  the  data  was  then  used  to  evaluate  the  flow  phenomenon  and  to  check 
the  validity  of  che  pretest  predictions  in  an  attempt  to  achieve  the 
objectives  of  thi/  program. 

5.1.2  Correlation  of  the  Data 

5. 1.2.1  Effect  of  Shock  Strength 

One  of  the  test  objectives  was  to  determine  the  effect  of 
shock  strength  on  the  loading  of  a cylinder  supported  horizontally  near 
the  ground.  Test  items  5.5a  and  5-5c  were  used  to  achieve  this  objec- 
tive and,  with  both  Shots  9 and  10,  yielded  four  sets  of  data.  These 
test  items  were  identical  in  size  (6-ft  diameter)  and  elevation  (5  ft 
above  the  ground).  The  values  of  the  shock  strengths  obtained  were: 

£ = 1.51  and  1.55  for  Shot  9>  and  ^ = 1.26  and  1.15  for  Shot  10  on 
items  5.5a  and  5 *50,  respectively.  These  values  of  shock  strength 
represent  a reasonable  spread  in  this  variable  for  loading  comparisons. 
(The  maximum  intended  shock  strength  was  approximately  ^ = 1.6.)  The 
basic  blast  wave  parameters  for  all  the  test  items  are  presented  in 
Table  4.2. 

The  pressure  data  were  plotted  in  normalized  form,  as  dis- 
cussed previously,  and  in  general  indicate  that  (a)  the  data  were  not 
sufficiently  accurate  to  determine  any  effects  of  shock  strength,  even 
during  the  diffraction  period,  and  (b)  the  pressure  magnitudes  in  many 
cases  appear  to  be  from  10  to  50  per  cent  low.  The  first  conclusion 
is  borne  out  by  examining  the  normalized  plots  presented  in  Figs.  5.3 
through  5*8.  It  was  not  possible  to  observe  any  consistent  or  inter- 
pretable trend  with  increasing  shock  strength  from  these  pressure-time 
plots.  This  randomness  of  behavior  prevailed  on  all  correlations 
which  were  made.  However,  in  only  a few  cases  did  the  plots  for  a given 
gage  position  differ  basically  in  shape. 

The  second  conclusion  follows  directly  from  the  discussion 
in  Section  5- 1-1  concerning  the  reliability  of  the  data.  In  view  of 
the  fact  that  the  strain  gage  data  provided  no  quantitative  results 
(see  Section  5*2),  the  test  objective  dealing  with  the  effect  of  shock 
strength  could  not  be  realized. 

Since  the  effect  of  shock  strer  -th  could  not  be  ascertained 
from  the  pressure  records,  the  data  were  averaged  and  further  reduced 
to  yield  what  is  believed  to  be  the  best  possible  estimate  of  the  load- 
ing on  a cylinder  one  radius  above  the  ground  (i.e.,  Ah/D  = 0.5). 

The  data  of  Figs.  5-5  through  5*8  were  averaged  linearly  in  most  cases; 
erratic  portions  of  the  records  were  ignored.  These  average  pressure- 
time records  were  also  modified  slightly  in  order  to  be  consistent 
with  corresponding  or  neighboring  records.  The  resultant  plots  are 
presented  in  Fig.  5*10  and  can  be  considered  to  represent  the  loading 
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on  a cylinder  one  radius  above  the  ground  and  struck  by  a shock  of 
strength  ^ = 1.52  (average  value  of  the  four  sets  of  data). 

The  average  pressure-time  plots  were  used  to  construct 
pressure  distributions  acting  on  the  cylinder  at  various  times  after 
shock  arrival.  This  procedure  was  carried  out  in  steps  of  one  r/U 
unit  up  to  a normalized  time  of  20  r/U  units  (approximately  50  ob). 

A few  typical  pressure  distributions  are  presented  in  Fig.  5»11* 

These  pressure  distribution  plots  were  then  used  to  compute  the 
average  pressure  in  the  horizontal  and  vertical  directions  (l.e.,  ratio 
of  total  force  to  projected  area).  The  net  horizontal  and  vertical 
loads  are  presented  in  Fig.  5.12.  The  magnitudes  found  are  question- 
able; however,  it  is  believed  that  the  time  details  of  the  loadings  are 
correct . 

The  significant  information  obtained  from  these  loadings 
are  the  clearing  and  buildup  times  for  the  respective  surfaces  (arcs). 
These  build  up  cmd  clearing  times  are  essentially  the  saiae  as  for  a 
two-dimensional  cylinder  in  free  space,  with  the  exception  of  the 
clearing  time  on  the  front  face.  The  predicted  loadings  for  a two- 
dimensional  cylinder  in  free  space  based  on  average  test  conditions 
is  given  in  Fig.  5- 15-  A plot  of  the  reflection  coefficient  referred 
to  in  Fig.  5-15  is  shown  in  Fig.  5*1^*  With  reference  to  the  experi- 
mental results  of  Fig.  5»10>  the  clearing  time  on  the  front  face  of 
a cylinder  located  one  radius  above  the  ground  is  estimated  to  be 
between  8 and  10  r/U  time  units.  As  seen  from  Fig.  5.15,  however, 
the  predicted  value  for  this  time  on  a cylinder  in  free  space  is  less 
than  5 r/U  units.  This  increase  is  reasonable,  however,  when  one  con- 
siders the  increase  in  pressure  due  to  ground  proximity  as  well  as  the 
greater  effective  clearing  distance  of  the  cylinder.  The  clearing  dis- 
tance of  the  cylinder  varies  from  the  order  of  r for  A h/D  ■ <» 

(i.e.,  in  free  space)  to  2r  for  A h/D  = 0 (l.e.,  touching  the  ground). 

5. 1.2. 2 Effect  of  Size 


Another  test  objective  deals  with  the  effect  of  size  of 
the  cylinder  on  the  loading.  Contributing  toward  this  objective  were 
Cylinders  5.3a  and  3.5b,  both  located  at  the  same  ground  range  and 
struck  by  identical  shock  waves  (^  ■ I.5I  for  Shot  9 and  £ a 1.26 
for  Shot  10).  These  cylinders  were  identical  in  geometry  but  differed 
in  absolute  size  by  a factor  of  4.  Cylinder  3.5a  was  6 ft  in  diameter 
and  located  5 ft  above  the  ground,  whereas  5.3b  was  I.5  ft  in  diameter 
and  located  9 in.  above  the  ground.  The  initial  Reynolds  numbers 
based  on  the  cylinder  diameters  are  shown  in  Table  5-1.  £/ 


2/  These  Reynolds  numbers  apply  to  flow  Just  behind  the  shock  front; 
the  flow  parameters  given  in  Table  4.2  have  been  computed  from 
the  conventional  Rankine-Hugonlot  theory  relating  the  state 
variables  across  the  shock  front. 


iSe 


a 

wsfmcfeB 


Table  5.1  - Initial  Reynolds  Numbers  for  Test  Cylinders 


Cylinder  5.3a 

R r,  M 14.2  X 10*^ 
eD 

for  Shot  9 

7.1  X 10^ 

for  Shot  10 

Cylinder  5.3b 

«eD  - 5-5  * 

for  Shot  9 

1.8  X 10^ 

for  Shot  10 

According  to  steady-state  concepts  (see  Aerodynamic  Drag  by  S.  F. 
Hoerner),  these  Reynolds  numbers  are  sufficiently  large  to  make  the 
drag  coefficient  relatively  independent  of  the  Reynolds  number  (critical 
Reynolds  number  S'  7 x 10^. 

The  simplified  pressure-time  plots  for  these  two  cylinders 
are  presented  in  Figs. 5. 15  through  5*20.  The  records  do  not  show  any 
large  deviations  in  the  drag  coefficient  due  to  a size  effect;  however, 
there  is  a considerable  and  generally  random  spread  in  the  values  of 
the  pressures.  This  spread  is  typical  of  all  the  pressure-time  records 
obtained  for  this  program.  The  actual  value  of  the  drag  coefficient 
could  not  be  determined  from  these  pressure-time  records. 

If  one  applies  a correction  factor  to  the  pressure  scales 
of  Fig.  5.15  such  that  the  peak  pressure  equals  1.82  p (o)  (i.e., 

the  predicted  value  for  a free  cylinder) , then  the  drag  coefficient 
on  the  front  face  is  approximately  -1,  and  the  drag  coefficient  on 
the  back  face  is  less  than  -h.  These  values  are,  of  course,  not 
realistic.  If,  on  the  other  hand,  one  assumes  that  the  pressure  scale 
of  Fig.  5.12a  is  correct,  then  the  drag  coefficients  are  even  more 
unrealistic. 

It  must  be  concluded,  therefore,  that  even  when  averaged, 
the  pressure  data  do  not  yield  any  information  about  the  actua."  value 
of  the  drag  coefficient.  While  it  was  realized  at  the  outset  that 
the  pressure  data  might  not  prove  sufficiently  accurate  to  determine 
drag  pressures  (see  Section  5.1) > this  is  nonetheless  unfortunate. 
Equally  unfortunate  is  the  fact  that  if  the  pressure  data  obtained  in 
this  test  are  as  reliable  as  one  can  expect  from  field  operations, 
the  determination  of  drag  pressures  associated  with  relatively  weak 
shocks  is  not  practical  by  means  of  direct  pressure  measurements  at 
the  present  time. 

The  obseivatlons  on  test  item  3.3b  tend  to  indicate  that  the 
occurrence  of  events  (i.e..  time)  scale  is  as  predicted.  The  response 
of  the  recording  system  was  Just  sufficient  to  detect  this  scaling. 

See  Figs.  5.15  through  5.20  for  these  trends.  The  magnitude  of  one 
r/U  unit  for  the  smaller  cylinder,  5.5b,  (r/U  Cf  0.6  ms)  was  too  small 
to  determine  clearly  the  exact  time  (±2r/U)  of  a particular  event. 
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5. 1.2. 3 Effect  of  Ground  Proximity 

A third,  and  perhaps  the  most  ijuportant,  objective  of  this 
test  was  to  determine  the  effect  of  ground  proximity.  The  test  items 
used  to  study  this  effect  were  3.3c,  3.5d,  and  3.3©.  These  cylinders 
were  identical  in  size  (6-ft  diameter)  and  were  elevated  3 ft,  1.5  ^t, 
and  4 in.  above  the  ground,  respectively,  resulting  in  Ah/D  values  of 
0.5,  0.25,  and  0.055.  The  test  items  were  located  at  the  same  ground 
range  and  were  subjected  to  essentially  identical  blast  waves  « 1.35 
for  Shot  9 and  ^ ■ 1.15  for  Shot  lO) • The  simplified  pressure-time 
records  are  presented  in  Figs.  5*21  through  5.26  for  these  cylinders. 

The  effect  of  ground  proximity  on  the  lo€«iing  of  a horizontal 
cylinder  was  not  determined  as  a result  of  this  test.  Determination 
could  not  be  made  for  two  primary  reasons.  First,  a majority  of  the 
pressure  records  (Gages  P4  and  P5,  Fig.  2.1)  which  would  contribute 
most  to  this  effect  were  either  lost  or  invalidated.  The  records  at 
gage  position  P5  (directly  facing  the  ground)  were  lost  for  test  items 
3.5d,  Shot  9>  and  5.5e,  Shots  9 and  10,  whereas  the  pressure-time  record 
on  test  item  3.5d,  Shot  10,  was  obviously  in  error  since  it  recorded 
a pressure  of  the  order  of  4p^(0)  for  a considerable  length  of  time 
(approximately  200  ms).  The  second  reason  is  that  these  particular 
gage  records  were  generally  poor,  both  with  respect  to  the  magnitude 
of  the  pressure  as  well  as  the  form  of  the  pressure-time  cuirves. 

(See  Figs.  5.21  through  5‘26.)  This  is  borne  out  if  comparisons  are 
made  of  the  records  from  Gages  PI,  P2,  P3,  and  (possibly)  P4  during 
the  initial  portion  of  the  loading  where  the  effect  of  the  reflecting 
surfaces  has  not  yet  influenced  the  pressure  at  these  gages.  Negative 
pressures  were  observed  initially  (see  Fig.  5-26)  at  the  pcfc..tlon  of 
Gage  P6  on  Cylinder  3.3e  during  both  Shots  9 and  10;  it  does  not  seem 
possible  that  the  initial  pressure  signal  could  actually  have  been 
negative.  There  are  many  corrections  or  processes  that  one  can  apply 
to  the  data  in  order  to  make  them  appear  more  realistic  (l.e.,  by  chang- 
ing scale  factors  or  by  shifting  time  scales,  etc.).  However,  in  all 
cases,  the  pressure  values  would  become  arbitrary  and  controversial, 
and,  it  is  felt,  this  procedure  would  not  Increase  the  meaningfulness 
of  the  test  results. 

5. 1.2. 4 Effect  of  Three-Dimensionality 

The  remaining  test  objective  concerned  three-dimensional 
effects  on  the  loading  (i.e.,  the  effect  of  finite  size).  This  objec- 
tive was  to  have  been  achieved  by  comijarlng  the  net  forces  obtained 
from  the  pressure  data  with  those  obtained  from  the  strain  gage  data. 

The  pressure  data  would  have  provided  two-dimensional  results,  at  least 
during  the  initial  loading  period,  whereas  the  strain  data  would  pro- 
vide the  three-dimensional  results.  As  discussed  in  later  sections, 
the  strain  data  did  not  yield  the  necessary  quantitative  Information 
to  effect  such  a comparison  and,  hence,  this  objective  could  not  be 
achieved.  This  is  indeed  unfortunate,  since  the  results  of  the  strain 
would  have  contributed  equally  as  well  to  all  test  objectives. 
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5.1.5  Reaults  of  Project  5.26.1  Pressure  Gages 


A portion  of  UPSHOT -KNOTHOLE  Project  5.26.1  (Tests  on  Effects 
of  POL  Installations,  WT-756),  is  related  to  the  present  program 
Inasmuch  as  both  deal  with  the  loading  on  structures  of  cylindrical 
shape.  There  was  one  vertical  cylindrical  storage  tank  in  Project 
5.26.1  (test  item  5.26bd)  which  was  instrumented  with  six  pressure 
gages  in  Shot  9*  (A  photograph  of  the  tank  is  shown  in  Fig.  2.5; 
the  pressure  gage  locations  are  Indicated  in  Fig.  2.6.  A description 
of  the  tank  is  given  in  Chapter  2.) 

The  tank  was  at  an  actual  ground  range  of  656O  ft  in  Shot  9 
and  was  struck  by  a clean  and  we 11 -developed  shock  wave.  The  peak 
overpressure  was  4.2  psi  and  the  positive  phase  duration  was  approxi- 
mately one  second.  All  of  the  pressure-time  records  exhibited  an 
appreciable  amount  of  "noise"  prior  to  shock  arrival,  and  it  must  be 
assumed  that  the  noise  level  did  not  change  thereafter.  Typical  linear- 
ized gage  records  (fast  playbacks)  are  shown  in  Fig.  5.27.  These 
pressure-time  observations  were  not  considered  in  the  final  report 
on  Project  5.26.1  and  are  discussed  in  this  section. 

The  roof  of  the  tank  collapsed  and  fell  into  the  tank  shortly 
after  shock  arrival.  This  time  has  been  estimated  to  be  about  100  ms 
from  motion  picture  films  taken  of  the  tank.  Thus,  the  observations 
were  invalidated  for  use  in  the  verification  of  loading  methods  on 
rigid  three-dimensional  cylinders  at  the  time  the  geometry  changed 
radically.  In  addition  to  the  change  in  geometry,  the  pressure  records 
from  the  side  gages  show  many  disturbances  which  are  not  considered 
to  be  due  to  flow  phenomenon  and  may  well  have  resulted  from  fragments 
of  the  roof  striking  the  gage  cables. 

In  general,  the  magnitudes  of  the  pressures  appear  to  be  low; 
however,  the  form  of  the  pressure-time  curves  for  the  gages  on  the  side 
of  the  tank  are  reasonable.  Pressure  Gage  P8,  which  was  facing  the 
blast  wave,  shows  a large  initial  peak  value  which  decays  rapidly  to 
a pseudo-steady-state  value.  Vfhile  the  magnitudes  were  probably  low, 
the  ratio  of  the  magnitudes  was  approximately  correct.  The  peak- 
reflected  pressure  cleared  in  approximately  15  ms,  which  is  nearly 
5r/U  time  units  and  agrees  with  the  clearing  time  for  a free  cylinder 
(see  Pig.  5.15).  This  gage  was  located  near  the  top  edge  of  the 
cylinder  and  indicates  that  the  three-dimensionality  of  this  structure 
does  not  affect  the  front  loading  appreciably. 

Pressure  Gage  PIO,  which  was  on  the  back  side  of  the  tank, 
indicated  a build  up  somevrtiat  more  rapid  than  predicted.  However,  this 
is  reasonable  since  the  gage  was  located  relatively  close  to  the  edge 
on  the  tank,  i.e.,  the  shock  wave  swept  around  from  the  side  in  addi- 
tion to  sweeping  over  the  top  edge.  This  record  also  shows  a vortex 
of  moderate  strength  (approximately  0.4p^(0))  sweeping  past  the  gage 
at  approximately  6r/U  time  units  after  the  shock  wave  had  reached 
the  back  of  the  tank.  The  observations  on  the  back  of  the  Project  5.5 
cylinders  indicate  a vortex  of  approximately  the  same  strength  sweep- 
ing past  at  from  4 to  5 r/U  time  units.  Pressure  Gage  P7,  which  was 
located  on  the  side  of  the  tank,  also  indicates  the  existence  of 
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a vortex  of  moderate  strength  (approximately  0.5p_(0)  ) at  the  time 
kr/U.  This  agrees  quite  well  with  the  observations  made  at  that  point 
on  the  Project  J.J  cylinders.  Pressure  Gage  P?  did  not  indicate  a small 
initial  peak  which,  it  is  believed,  should  have  existed. 

The  records  of  roof  Gages  PI,  PJ,  and  P6  do  not  seem  to  be 
realistic,  and  since  the  roof  failed  early,  the  records  do  not  provide 
useful  information, even  in  the  very  early  stages  prior  to  roof  failure. 

5.2  STRAIN  MBASUREMEMTS 


5.2.1  Introduction 


The  analysis  and  discussion  of  the  strain  data  is  presented 
in  the  following  sections.  Section  5.2.2  is  concerned  with  a dis- 
cussion of  the  frequency  response  of  the  cylinder- sensor  system,  and 
an  attempt  is  made  to  identify  the  physical  origin  of  the  various  fre- 
quencies observed  in  the  records.  It  appears  likely  that  a well-defined 
low-frequency  oscillation  found  in  practically  all  records  may  be  due 
to  a rigid-body  motion  of  the  cylinder  foundation  In  the  soil.  In  any 
event,  the  frequency  response  of  the  test  items  is  considerably  more 
complicated  than  originally  anticipated,  and  it  is  not  considered 
feasible  to  determine  the  response  of  this  system  to  even  a known  input 
by  analytical  means. 

A comparison  of  peak  horizontal  force  as  computed  from  the  strain 
data  and  pressure  data,  together  with  the  predicted  values  for  an  iso- 
lated cylinder,  is  given  in  Section  5 *2. 5-  Based  on  this  comparison, 
the  strain  data  yield  results  which  are  in  all  cases  lower  than  expected, 
and  in  some  cases  lower  than  is  reasonable.  The  vertical  force  acting 
at  the  time  of  maximum  horizontal  force  indicates  a random  directional 
behavior.  That  is,  in  some  cases  the  vertical  force  appears  to  act  up- 
wards; and  in  others,  downwards.  The  fact  that  an  initial  downward 
force  is  considered  unrealistic  points  to  either  an  uncertainty  in  the 
measured  peak  strains  of  about  30  per  cent  or  more,  or  the  fact  that 
this  interpretation  of  the  records  is  invalid. 

The  magnitude  and  time-details  of  the  incident  blast  loading 
could  not  be  adequately  determined  from  visual  inspection  of  the  strain 
records  alone.  That  this  might  be  the  case  was  anticipated  from  the 
outset,  and  a numerical  method  of  analysis  designed  to  determine  the 
input  (i.e.,  the  blast  loading)  from  the  known  output  (i.e.,  the  strain 
data)  was  devised  during  the  pretest  planning  phase.  This  method  of 
data  reduction,  referred  to  as  transient  analysis,  was  applied  to  all 
strain  data  which  appeared  reasonable,  but  for  one  reason  or  another, 
the  results  of  this  analysis  proved  largely  unsatisfactory.  The  develop- 
ment of  the  method  of  transient  analysis,  together  with  a critical 
discussion  of  its  general  applicability  as  it  is  currently  understood, 
is  presented  in  Appendix  A;  the  actual  work  carried  out  and  the  results 
obtained  are  given  in  Section  5«2.U. 
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A first  application  of  the  method  of  transient  analysis  proved 
inconclusive.  A second  application  was  begun  but,  since  it  did  not 

seem  as  if  the  necessary  additional  work  required  to  prove  or  disprove  .j 

the  method  in  the  present  application  would  be  consistent  with  the  ' 

scope  of  the  program,  this  work  was  not  completed.  Thus,  no  definite  ■] 

conclusions  have  been  reached  as  to  the  feasibility  of  this  method  in  jj 

general,  although  it  must  be  admitted  that  much  has  been  learned  regard-  |i 

Ing  certain  aspects  of  the  analysis  and  the  general  complexity  of  the  j 

computational  procedures  \rtiich  was  not  at  all  appreciated  at  the  outset.  |j 

And  it  must  further  be  admitted  and  recognized  that  the  problem  of 
data-reductlon  stemming  from  so-called  net  force  measurement  systems 
such  as  in  the  present  case  (l.e.,  where  the  test  items  represent  a 
fairly  complicated  dynamic  system)  is  still  largely  an  unsettled  issue. 

In  view  of  the  unsuccessful  results  of  the  numerical  data- 
reduction  scheme,  it  was  necessary  to  follow  another  approach  in  an 
effort  to  obtain  the  desired  information  from  the  strain  data.  The 
method  adopted,  which  can  be  called  only  a qualitative  approach  to  the 
problem,  was  an  attempt  to  remove  the  harmonic  ccxnponents  of  the  res- 
ponse by  lnsi)ection  and  thereby  yield  at  least  the  shape  of  the  forcing 
function.  A discussion  of  this  work  and  the  rather  conflicting  results 
obtained  is  given  in  Section 

In  connection  with  the  qualitative  approach  indicated  above, 
a brief  effort  was  directed  toward  utilization  of  the  ARF  analog  com- 
puter as  a curve-fitting  device.  The  work  is  described  in  Section 
5.2.6.  While  no  conclusive  results  were  obtained,  the  technique  appears 
to  be  a promising  one. 

5.2.2  Frequency  Responee  of  Cylinders 

A typical  strain-time  record  is  shown  in  Fig.  5>28a.  Even  a 
brief  inspection  of  this  record  is  sufficient  to  point  out  that  the 
cylinder-sensor  arrangement  did  not  respond  as  what  could  be  termed 
a simple  dynamic  system.  Quite  the  contrary,  as  at  least  four  princi- 
pal frequencies  can  be  identified  in  the  majority  of  the  records. 

An  attempt  is  made  in  this  section  to  identify  these  various  frequencies 
so  as  to  Interpo^t  better  the  strain  data. 

Average  values  of  the  principal  periods  as  determined  from  s 
visual  inspection  of  the  strain  records  are  shown  in  Table  5*2.  Where 
no  periods  are  shown,  either  no  usable  record  was  obtained  (see  Table  U.3), 
or  the  periods  were  not  clearly  defined.  While  the  results  of  Table  5-2 
could  no  doubt  be  improved  by  a statistical  approach  to  the  data  (the 
autocorrelation  function  was  calculated  for  one  set  of  data),  these 
approximate  values  are  adequate  for  the  present  discussion.  It  might 
also  be  mentioned  that  a much  lower  frequency  of  about  5 cps  (a  period 
of  200  ms  or  greater)  not  shown  in  Table  5*2  was  evident  in  a number 
of  the  original  record  playbacks.  The  discrepancies  noted  in  frequencies 
for  the  same  cylinder  in  the  two  shots  may  either  be  due  to  the  inade- 
quacy of  determining  these  values  by  Inspection,  or  to  the  differences 
in  excitation  brought  about  by  the  difference  in  orientation  of  the 
shock . 
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The  strain  records  are  direct  evidence  of  the  highly  involved 
response  of  the  cylinder-sensor  system,  which  is  that  of  a relatively 
thin  circular  shell  with  restrained  ends  on  flexible  supports.  A de- 
tailed dynamic  analysis  of  this  system  was  not  carried  out  in  view  of 
the  technical  difficulties  and  uncertainties  involved.  In  fact,  it  is 
precisely  one's  lack  of  ability  to  characterize  adequately  the  response 
of  the  test  setup  that  makes  the  analysis  and  interpretation  of  the 
strain  records  such  a formidable  problem.  For  the  purpose  of  identify- 
ing the  principal  frequencies,  however,  it  is  possible  to  consider 
certain  idealized  response  modes  of  the  cylinder-sensor  system  indi- 
vidually. Several  such  idealized  response  modes  and  associated  fre- 
quencies (and  periods)  are  indicated  in  Table  5*3 • The  entries  in 
this  table  are  discussed  in  the  following  paragraphs. 

If  the  cylinder  were  sufficiently  rigid,  it  is  reasonable  to 
expect  that  the  principal  symmetrical  response  mode  of  the  system  would 
be  translation  of  the  cylinder  parallel  to  the  plane  of  the  sensor 
supports.  If  the  bending  resistance  of  the  sensor  bars  is  neglected 
as  compared  to  their  axial  stiffness,  this  motion  can  be  idealized  as 
that  corresponding  to  a single-degree-of-freedom  mass-spring  system. 

This  is  the  first  mode  listed  in  Table  5*3>  the  required  data  are 
given  in  Chapter  2,  Table  2.1.  Two  extreme  values  for  the  period  of  this 
mode  are  given  for  both  the  6-  and  1.5-ft  diameter  cylinders.  This 
variation  is  due  to  the  uncertainty  regarding  the  effective  length  of 
the  sensor  bars  in  this  motion  (see  note  in  Table  2.1).  The  higher 
frequency  corresponds  to  the  unsupported  length  of  the  bars;  the  lower 
frequency  is  based  on  the  over-all  length  of  the  bars  (see  Fig.  2.2). 

The  second  mode  considered  is  that  of  simple  beam  vibration  of 
the  cylinder  with  the  sensors  acting  as  rigid  supports,  but  providing 
no  end  moment  restraint.^  This  frequency  is  computed  frcji  the  well- 
known  relation  for  a simply  supported  beam  and  Incorporates  first-order 
corrections  for  shear  deformation  (see  Vibration  Problems  in  Engineering 
by  S.  Timoshenko).  The  frequencies  associated  with  this  mode  are  about 
twice  as  great  as  the  rigid  cylinder  motion. 

Another  idealized  mode  is  that  of  the  vibration  of  an  isolated 
thin  circular  shell.  In  the  present  case  this  is  a rather  complicated 
problem,  even  if  the  influence  of  the  sensor  supixjrts  is  neglecced, 
since  the  effect  of  the  end  plates  is  to  maintain  a fixed  diameter. 
Solutions  for  the  vibration  of  a thin  cylindrical  shell  of  finite  length 
whose  ends  remain  circular  but  otherwise  unrestrained  are  available 
(see  Mathematical  Theory  of  Elasticity  by  A.  E.  H.  Love).  The  frequencies 

5/  The  frequencies  of  vibration  of  a beam  mounted  on  the  same  type 

of  sensor  supports  and  including  the  flexibility  and  moment  restraint 
of  the  supports  was  carried  out  in  connection  with  a pretest  analy- 
sis on  Operation  TEAPOT  Air  Force  Project  5*2  (Study  of  Drag  Loading 
of  Structures  in  and  out  of  the  Precursor  Zone,  WT-1124).  The 
results  obtained  did  not  differ  too  greatly  from  the  frequencies  of 
a conventional  simply  supported  beeun. 
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TABLE  5 >3  - Suranary  of  Frequency  Response  of  Cylinders 


Assumed  Mode  of  Response 

Frequency  f : Period  T 

Cylinders 

3.5a,  c,  d,  e 

Cylinder 

5.3b 

1.  Rigid  Cylinder  on  Flexible 
Supports 

67  < f < 99  cps 

10  < T < 15  mspc* 

155  < f < 199  cps 

5 < T < 7 mspc 

2.  Cylinder  as  Simple  Beam 
on  Rigid  Supports 

f = 122  cps 

T ■ 8 mspc 

f = 480  cps 

T = 2 mspc 

5.  Cylinder  as  Thin 

Circular  Shell 

(a)  Extensional  Modes 

(b)  Inextensional  Modes 

1st  Mode 

Ifth  Mode 

f = 7^6  cps 

T ■ 1.4  mspc 

f « 10  cps 

T ■ 1000  mspc 

f » 85  cps 

T = 12  mspc 

f = 154  cps 

T = 6. 5 mspc 

1^.  Rigid  Body  Motion 
in  Soil 

(a)  Rotational  Mode 

(b)  Vertical  Translational 

Mode 

(c)  Horizontal  Transla- 

tional Mode 

f » 16  cps 

T ■ 62  mspc 
f ■ 12.8  cps 

T ■ 78  mspc 
f = 6.5  cps 

T » 155  msps 

raspc  = milliseconds  per  cycle 
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5.2.5  Comparison  of  Peak  Transmitted  Forces 


It  is  of  considerable  interest  to  determine  the  peak  forces  act- 
ing on  the  test  cylinder  as  Indicated  by  the  strain  records.  The  initial 
peak  strains  were  read  directly  from  the  linearized  strain  records  and 
converted  to  total  force  per  unit  projected  area  of  the  cylinder  trans- 
mitted to  the  sensor  bars  according  to  Equation  5 .3-  These  results, 
in  the  form  of  horizontal  and  vertical  force  components,  along  with 
the  peak  strain  date^  are  shown  in  Table  5*^.  It  might  be  noted  that 
the  strain  gage  system  was  apparently  wired  so  that  a tensile  force 
in  the  sensor  bars  was  recorded  as  a negative  strain.  These  signs 
have  been  reversed  in  Table  5.^  so  that  a positive  F.  indicates  a 
horizontal  force  directed  downstream,  and  a positive  r Indicates 
a vertical  force  directed  upward. 

According  to  Equation  5.3  the  strains  in  the  front  and  rear 
sensor  bars  should  be  numerically  equal  (but  of  opposite  sign)  in  the 
absence  of  a vertical  force  component  for  the  static  case.  The  majority 
of  the  strain  recosrds  do  indeed  show  a well-defined  initial  peak  strain 
which  is  of  the  same  order  of  magnitude  for  both  front  and  rear  sensors 
of  the  same  cylinder.  Thus,  in  those  cases  where  the  data  from  only 
one  sensor  bar  were  obtained,  the  peak  horizontal  force  indicated  in 
Table  5.4  is  based  on  the  assumption  that  the  same  magnitude  strain 
occurred  in  the  other  sensor  bar.  Due  to  the  uncertainty  of  this 
assumption,  however,  the  vertical  force  component  which  depends  on  the 
difference  of  the  magnitudes  of  peak  strains  was  not  computed  unless 
records  fran  both  sensor  bars  were  obtained. 

The  measured  reactions  (i.e.,  F.  and  F ) in  general  cannot  be 
compared  directly  with  the  applied  force “due  to  tSe  influence  of  the 
dynamic  characteristics  of  the  cylinder-sensor  system  on  these  reactions. 
In  the  present  case,  however,  the  appearance  of  well-defined  initial 
peaks  suggests  the  following  approximate  scheme  for  deducing  the  peak 
applied  force  from  the  measured  peak  reactions.  It  is  assumed  that 
during  the  first  half-cycle  of  motion  (i.e.,  in  the  time  of  rise  to  the 
initial  peak)  the  cylinder  acts  as  a single-degree-of-freedom  system 
with  period  equal  to  twice  the  observed  rise  time.  Knowing  this  period 
and  assuming  undamped  motion  initially,  it  is  possible  to  r'late  the 
dynamic  overshoot  of  the  system  to  the  peak  applied  loading  jy  means  of 
the  so-called  dynamic  load  factor  (DLF)^  for  the  appropriate  type  of 
forcing  function. 

The  average  rise  time  for  the  6-ft  diameter  cylinders  was  found 
from  the  strain  records  to  be  about  8 ms,  indicating  an  initial  period 
of  about  16  mspc;  for  the  quarter-scale  cylinder,  this  period  was  found 
to  be  about  6 mspc  (note  fresn  Table  5*2  that  such  periods  were  found 
throughout  the  entire  history  of  the  motion) . According  to  the  pretest 
load  prediction  scheme  for  a cylinder  in  free  space  (see  Fig.  5.13), 
the  net  horizontal  load  during  the  diffraction  phrase  rises  linearly 
from  zero  to  a maximum  in  one  r/U  unit  (about  2.4  and  0.6  ms  for  the 
6-  and  1.5-ft  diameter  cylinders,  respectively)  and  then  decays  linearly 

4/  Peak  dynamic  response  = DLF  x (static  response  to  peak  applied  load). 
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to  pseudo-steady-state  In  between  4 and  5r/U  units.  Inas.-nuch  as  this 
decay  time  is  somewhat  less  than  the  initial  periods,  it  is  sufficient 
to  consider  only  this  triangular  diffraction  loading  in  estimating  the 
desired  DLF.  From  solutions  available  in  the  literature  (Elastic 
Structures  Under  Rapidly  Applied  Loading  - III,  Mechanical  World  and 
Engineering  Record,  No.  257^),  it  is  found  that  DLF  =1.5  for  the 
6-ft  diameter  cylinder  and  DLF  = 1.5  for  the  1.5-ft  cylinder  for 
both  shots.  Accordingly,  the  column  labeled  Fjj/DLF  in  Table  5-4 
represents  the  estimated  peak  applied  horizontal  force  as  indicated 
by  the  strain  records. 

The  measured  values  of  Pj^  shown  in  Table  5*4  are  taken  from 
the  results  of  Fig.  5.12a  for  the  cylinder  configuration  Ah/D  = 0.5, 
where  the  peak  pressure  was  found  to  be  about  1.3Pq(o).  The  predicted 
values  of  F^j  shown  in  the  table  were  computed  from  the  relation 

F 

[c=yl  hor  - l]  * i (M) 

which  is  applicable  to  a cylinder  in  free  space  (see  Fig.  5.15).  A plot 
of  the  average  reflection  coefficient,  C hor’  given  in  Fig.  5-14. 

With  reference  to  Table  5.4,  it  is  seen  that  the  peak  force  as 
determined  by  the  strain  gage  data,  in  every  case,  appears  to  be  low; 
in  most  cases,  even  less  than  the  side-on  pressure  acting  at  that  time. 
The  values  of  F (i.e.,  the  vertical  force  component  at  the  time  the 
peak  Fj^  occurs)  also  appear  to  be  unrealistic.  More  often  than  not, 
the  results  indicate  a net  vertical  force  acting  downward  on  the  cylin- 
der, which  is  contrary  to  all  previous  experience.  And  there  is  no 
consistency  in  the  apparent  direction  of  Fy  between  the  results  of 
Shots  9 and  10.  Although  the  magnitude  of  Fy  increases  with  increas- 
ing A h,  as  would  be  expected,  this  trend  is  not  consistent  between 
the  two  shots. 

In  view  of  the  low  values  of  Fjj,  as  computed  from  the  strain 
data,  and  the  randomness  in  sign  of  Fy,  it  appears  that  either  the 
association  between  the  measured  strains  and  the  applied  forces  is 
invalid,  or  that  the  magnitudes  of  the  strain  data  are  unreliable. 

It  is  impossible  to  assess  fully  the  validity  of  converting  the  peak 
measured  strains  to  applied  forces  short  of  undertaking  a comprehensive 
dynamic  analysis  of  the  system.  As  long  as  the  cylinder-sensor  system 
responds  initially  in  only  one  mode,  the  use  of  a DLF  is  valid.  While 
the  proper  value  of  the  DLF  is  open  to  some  question,  even  a lower  bound 
value  of  DLF  =15/  results  in  peak  forces  which  are  still  considered 
low.  However,  it  is  possible  that  certain  harmonic  components  of  the 
dynamic  system  act  so  as  to  appreciably  distort  the  initial  time  details 
of  the  loading.  In  such  case,  the  initial  peak  strain  might  no  longer 
be  proportional  to  the  applied  load,  and  the  comparison  of  Table  5*4 
would  be  invalid.  At  later  times  in  the  response,  however,  these  com- 
ponents may  be  damped  out,  and  the  shape  and  magnitude  of  the  applied 
force  might  be  readily  obtained. 

^ DLF  =1  is  a lower  bound  value  for  the  periods  and  duration 
of  diffraction  loading  considered  here. 
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The  strain  data  for  Cylinder  3 •3a,  Shot  10,  tend  to  support 
this  view.  As  indicated  in  Table  5.h,  the  peak  recorded  force  of 
0.45  psi  is  felt  to  be  low  by  a factor  of  10  or  more.  However,  as 
discussed  in  Section  5*2.4,  after  about  20  ms  the  strain  data  appear 
quite  realistic  and,  in  fact,  supports  steady-state  wind-tunnel  results 
for  drag  coefficients.  Similarly,  the  randomness  in  behavior  of  ' ost 
of  the  other  cylinders  could  be  explained  in  terms  of  a complicated 
response.  In  Shot  9,  identical  Cylinders  3 *36  and  c were  in  regions 
of  different  shock  strength,  yet  the  measured  peak  strains  were 
practically  identical.  Cylinders  3*3c,  d and  e had  decreasing  Ah/D 
in  the  order  given  (see  Fig.  2.1)  and  were  exposed  to  the  same  loading. 
However,  the  peak  horizontal  force  does  not  show  a trend  consistent 
with  the  variation  in  Ah/D,  in  that  the  intermediate  cylinder,  3.3<1, 
in  both  Shots  9 and  10  indicated  the  higher  loading.  From  the  other 
point  of  view,  however,  the  results  for  Cylinders  3 *50  and  d in  Shot  9, 
and  3*3^  and  e in  Shot  10  show  extremely  well-defined  initial  peaks 
(see  Fig.  5«28)  and  it  seems  reasonable  that  the  present  interpretation 
is  valid,  at  least  in  these  cases. 

As  seen  in  Table  5-4,  the  peak  strains  in  the  two  sensor  bars 
differ  in  one  case  by  as  little  as  11^  in. /in.  and  in  another  by 
35/xln./in.  Under  static  load  conditions  the  reliability  of  the  SR-4 
strain  gages  used  in  the  test  is  only  about  ±20^ in. /in.  Thus  the 
inconsistency  in  the  vertical  force  could  be  explained  in  part  by 
reasonable  inaccuracies  in  the  strain  data.  Of  course,  it  is  possible 
that  errors  in  strain  or  strain  calibration  could  account  for  all  the 
uncertainty  noted  here.  It  would  tjxen  follow  from  a consideration  of 
F^  that  the  uncertainty  in  strain  magnitude  is  of  the  order  of 
30  per  cent  — or  an  uncertainty  in  pressure  of  from  1 to  2 psi. 

However,  it  should  be  noted  that  the  magnitude  of  the  vertical  force 
shown  in  Table  5*4  in  some  cases  is  not  inconsistent  with  what  was 
found  from  the  pressure  gage  data.  For  example  the  results  of  Fig.  5 •12b 
(applicable  to  Cylinder  3 *30)  indicate  a net  upward  force  of  about 
0.7  pai  for  that  cylinder  (in  Shot  9)  at  the  time  the  peak  horizontal 
force  occurred.  Considering  the  fact  that  there  will  be  some  time  lag 
in  the  response  of  the  cylinder  to  this  load,  a measured  value  of 
0.34  psi  is  not  unrealistic. 

There  is  an  alternative  explanation  to  the  apparent  error  in 
sign  of  Fy  and  low  values  of  the  raaocimum  horizontal  force.  This  is 
that  the  recording  equipment  may  have  missed  the  initial  peak  strain 
in  some  cases,  and  what  has  been  accepted  as  the  first  maximum  of  the 
strain  is  in  reality  the  first  minimum  which  occurs  one-half  cycle 
later.  Hot  only  would  this  account  for  the  low  values  of  Fjj,  but  it 
would  imply  a change  in  sign  of  F^,  since  what  was  assumed  to  be  posi- 
tive strain  now  takes  on  a negative  value.  (The  signs  would  then  agree 
with  BRL's  sign  convention  on  Che  original  playbacks.)  However,  in 
view  of  the  rated  response  time  of  the  recording  equipment  as  well  as 
in  inspection  of  the  original  record  playback,  this  explanation  does 
not  appear  to  be  realistic. 
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To  sui.'nr.arize,  the  *naxiiiiurn  recorded  stratus  point  to  a general 
randc:''ness  and  inconsistency  in  behavior  when  the  results  are  taken 
as  a whole.  Whether  the  interpretation  of  these  data  is  invalid  due  to 
the  complexity  of  the  dynamic  system,  or  whether  there  exists  a large 
uncertainty  In  the  records  themselves  due  to  errors  in  'iiafn Ltude  or 
calibration,  cannot  definitely  be  stated;  some  evidence  is  presented 
for  either  reason. 

In  retrospect  it  is  now  evident  that  in  field  tests  of  this 
nature,  as  much  information  as  possible  should  be  obtained  regarding 
the  dynamic  systems  being  tested.  A minimum  of  such  information  would 
consist  of  pretest  and  post-test  static  pull  tests  which  would  serve 
to  define  the  linear  range  of  the  system  and  check  out  the  strain  gage 
system  and  calibration.  Some  sort  of  dynamic  tests  should  also  be 
considered  in  order  to  determine  the  principal  frequencies  and  damping 
of  the  cylinder-sensor  system.  Information  of  this  type  would  undoubt- 
edly have  aided  in  the  study  and  interpretation  of  the  strain  data. 

5.2.4  Graphical  Analysis  of  the  Strain  Records 

The  results  of  a rather  comprehensive  visual  inspection  of  the 
individual  strain  records  is  presented  in  this  section.  For  lack  of 
a better  name,  the  approach  followed  is  referred  to  as  graphical  analy- 
sis, and  it  should  be  emphasized  at  the  outset  that  the  method  and 
results  obtained  are  only  .meaningful  in  a qualitative  sense.  Briefly, 
the  method  consisted  of  combining  the  two  strain  records  for  each 
cylinder  into  horizontal  and  vertical  components  (where  possible),  and 
constructing  what  appeared  to  be  the  most  reasonable  median  curve  about 
which  the  combined  records  appeared  to  be  oscillating.  This  median 
curve  was  obtained  by  first  estimating  the  envelopes  to  the  record  and 
then  constructing  a curve  which  passed  midway  between  the  upper  and 
lower  bounds.  This  procedure  was  repeated  several  times  until  a reason- 
ably smooth,  nonoscillating  curve  was  obtained;  this  curve  was  then 
taken  as  representing  the  general  shape  of  the  forcing  function  being 
considered.  In  effect,  the  graphical  approach  filters  out  the  harmonic 
components  of  the  response  and,  in  principle,  accomplishes  what  was 
intended  with  the  method  of  transie.nt  analysis. 

Before  considering  the  results  obtained,  there  are  several 
points  to  be  mentioned  in  connection  with  the  method  of  analysis  and 
the  records  themselves.  The  method  of  transient  analysis  (discussed 
in  the  following  section  and  in  Appendix  A)  was  applied  to  individual 
strain  records,  since  a preliminary  inspection  of  the  data  indicated 
some  uncertainty  in  the  strain  calibration  and  in  the  locat  Lon  of  shock 
arrival  time.  It  was  decided  to  base  the  present  approach,  however, 
on  the  combined  data,  (i.e.,  horizontal  and  vertical  strain  components) 
inasmuch  as  a direct  comparison  with  the  results  of  the  pressure  gage 
analysis  and  predicted  values  was  desired. 

As  discussed  in  connection  with  the  strain  data  for  Cylinder 
5.5a  on  Shot  10,  the  results  obtained  do  not  appear  to  be  seriously 
affected  within  the  range  of  uncertainty  in  locating  a consistent 
zero  time  between  the  two  strain  records  for  a particular  cylinder. 
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In  those  cases  where  well-defined  initial  peak  strains  occurred, 
a common  zero  time  was  established  by  matching  these  peaks  on  the  two 
records.  In  those  cases  where  the  initial  peaks  were  not  well-defined, 
the  zero  time  was  determined  in  a somewhat  arbitrary  fashion.  In  any 
event,  the  results  are  not  believed  to  be  materially  affected  by  errors 
in  the  location  of  zero  time.  With  respect  to  possible  errors  or 
Inconsistencies  in  the  strain  calibration,  little  can  be  said.  However, 
as  mentioned  in  connection  with  a comparison  of  peak  force  values, 
the  general  inconsistency  of  the  results  obtained  may  in  part  reflect 
on  the  validity  of  the  strain  data  itself. 

As  indicated  in  Table  4.5,  the  majority  of  the  linearized 
strain  plots  extended  to  only  about  500  ms  after  shock  arrival.  Since 
the  duration  of  the  blast  was  about  1 sec  for  all  cylinders,  the 
linearized  data  provide  information  for  only  about  the  first  third  of 
the  loading  period.  However,  in  view  of  the  nature  of  the  results 
obtained,  this  is  not  felt  to  be  a serious  limitation,  inasmuch  as 
the  (nonlinear)  slow  playbacks  provide  adequate  information  in  those 
cases  where  details  at  later  times  are  of  interest. 

Normalized  overpressure  and  dynamic  pressure-time  cuirves  for 
Shots  9 and  10  are  shown  in  Fig.  5*29*  The  overpressure  curves  are 
based  on  the  commonly  accepted  empirical  relationship: 

P^(t)  = pjo)  e ° (l-t/t^)  (5-2) 

where  P^(0)  is  the  incident  overpressure,  t^  is  the  duration  of  the 
positive  phase  of  the  blast,  and  c is  an  empirically  determined  constant. 
In  general,  c is  known  to  depend  on  both  scaled  height-of -burst  and 
ground  range  (UPSHOT -KNOTHOLE  Project  5'1,  Tests  on  the  Loading  of 

Building  and  Equipment  Shapes,  WT-721).  However,  for  the  present  dis- 
cu8sion7  c = ^ and  c = 1 are  found  to  represent  adequately  the  over- 
pressure wave  form  at  both  test  locations  in  Shots  9 and  10,  respectively. 

While  the  relationship  between  dynamic  pressure  and  overpressure 
at  the  shock  front  is  well-known,  the  wave  form  of  the  dynamic  pressure 
curve  is  considerably  less  certain,  and  adequate  measurements  are  not 
available  at  the  test  locations.  In  the  past,  it  has  been  commonly 
assumed  that  the  Rankine-Hugoniot  conditions  applicable  at  the  shock 
front,  apply  throughout  the  wave.  For  small  shock  strengths,  as  in 
the  present  case,  this  is  tantamount  to  taking  the  dynamic  pressure 
time-cuirve  proportional  to  the  square  of  the  overpressure  time-curve. 
While  the  validity  of  this  approach  may  be  questioned,  the  d’^namic 
pressure  curves  shown  in  Fig.  5-29  are  based  on  such  an  assumption. 
Inasmuch  as  the  application  of  these  curves  is  for  a qualitative  com- 
parison with  the  results  of  the  analysis  discussed  in  this  section, 
the  present  approach  is  considered  acceptable. 

The  results  obtained  by  the  so-called  method  of  graphical 
analysis  for  the  six  sets  of  strain  data  considered  are  summarized  in 
Table  5*5>  which  briefly  indicates  the  character  of  the  forcing  func- 
tions so  determined. 
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T^BLB  5*5  - Sumoaxy  of  Results  Obtained  by  Graphical  Analysis 


General  Shape  of  Forcing  Function 

Cylinder 

Horizontal  Conponent, 

Vertical  Component,  Fy 

j.5c.  Shot  9 

(See  Flge,  5-16,  5.J7) 

Ah  - 1 

p^O)  ■ 4.6  pel 

Force  directed  upetrean  after 
about  30  ne.  At  28o  oa  force 
le  about  0.6  pel;  dynamic 
preaaure  at  thle  time  la 
about  0.08  pel.  Apparently 
a eerloua  baaeline  shift. 

Force  directed  upward  which 
tends  to  about  0,6  pal  at  about 

280  as.  If  baseline  shift  of 
this  amount  were  Included, 

Fy  would  appear  to  oscillate 
about  zero. 

}.)d,  Shot  9 

(See  Flge.  5.}8,  5.39) 

Ah  • ^ 

p^O)  » 4.6  pel 

Shape  of  Fjj  uncertain  between 

20  to  100  me.  Indicates  ceak 
unit  force  of  about  one  psl 
directed  upstream  at  60  ne. 
Between  100  and  200  as,  Fj^ 
oscillates  about  curve  having 
approximate  wave  fom  as  dy- 
namic pressure  curve;  "drag 
coefficient"  about  1.3  during 
this  tine. 

Force  oscillating  about  zero 
or  possibly  a slightly  negative 
(l.e.,  upMird)  value. 

3.3«,  Shot  10 
(See  Flge.  5.42, 

5.1*3,  5.44) 

Ah  • j 

Pg](0)  • 5.5  pel 

No  appearance  of  Initial  peak 
force.  (See  Table  5.1i>) 

Indicates  wave  fora  practi- 
cally parallel  to  dynamic 
pressure  curve  fron  about  15  to 
150  ms.  Average  'drag 
coefficient"  is  about  0.27. 

Upward  force  gradually  Increasing 
from  zero  to  about  0,15  psl  at 

40  as  and  tends  to  about  0.21  psl 

In  140  as.  Dynamic  pressure  about 
0.17  psl  at  latter  time.  Rate  of 
Initial  rise  may  be  more  rapid 
depending  on  choice  of  zero  time. 

J.Jb,  Shot  10 

Ah  • § 

PgiO)  m 5.5  pel 

Inconclusive;  no  Indication 
of  sustained  Initial  peak. 

F|j  appears  to  oscillate  about 
tero  or  possibly  a slightly 
negative  (l.s.,  upstream) 
value. 

Fy  very  similar  to  that  of  3.3*, 
Shot  10;  rises  to  about  0.I6  psl 

In  30  ms  and  tends  to  about  0.23 
psl  In  160  ms.  Dynamic  pressure 

Is  about  0.15  psl  at  latter 
time. 

5.5d,  Shot  10 

(See  Flge.  5.40,  5.1*1) 

Ah  - ^ 

PglO)  • 2.0  pel 

Fh  very  similar  In  shape  to 
that  obtained  for  same  cylin- 
der In  Shot  9.  Shape  uncertain 
between  about  30  and  100  ms. 
Between  about  140  and  240  ms, 

P{,  sopears  to  be  well  defined, 
Indicating  a curve  which  decays 
somewhat  more  rapidly  than  does 
the  dynamic  pressure  curve  dur- 
ing this  time;  the  "drag  coef- 
ficient" varies  from  about  5.O 
at  160  ms  to  1.4  at  230  ms. 

Indicates  force  oscillating 
about  zero. 

3.3ei  Shot  10 

(See  Flge.  5-45,  S.Ue) 

Ah  ■ ^ • 0.33 

Pglo)  m 2.0  pel 

Fjj  similar  to  3.3d,  Shot  10, 
although  decay  of  curve  between 
about  100  and  200  ms  not  as 
well  defined.  The  sverage 
"drag  coefficient"  appears  to 
be  about  6. 

Initial  rise  not  well  defined. 
Decays  to  about  0.15  psl  at  I60  as 
along  a curve  reasonably  parallel 
to  dynamic  pressure  curve,  between 
about  Bo  to  160  ms;  magnitude  of 

Fy  Is  from  2.5  to  3 times  the 
dynamic  pressure  during  this  time. 

Referring  to  Table  Fig.  5*50,  the  results  for  Cylinder 

5.5c,  Shot  9>  indicate  a strong  possibility  that  an  appreciable  base- 
line shift  occurred  in  the  strain  record.  The  original  playback  of  one 
strain  gage  (5.3cSl)  suggests  a baseline  shift  at  about  200  ms  after 
shock  arrival;  the  other  record  shows  no  peculiarities  at  this  time. 
However,  since  both  slow  playbacks  end  at  about  600  ms  after  shock 
arrival  (i.e.,  prior  to  the  duration  of  the  blast)  no  information  is 
available  regarding  the  condition  of  the  baseline  at  later  times. 

In  view  of  the  fact  that  Fj^  as  determined  from  the  data  appears  to 
be  the  only  reasonable  interpretation,  especially  at  later  times,  the 
data  for  this  cylinder  are  considered  to  be  erroneous.  Vihlle  the  res- 
ponse during  the  initial  loading  period  may  well  be  meaningful,  it  is 
not  possible  to  deduce  the  necessary  detail  of  the  loading  function 
by  the  present  method. 

Referring  to  Table  5-5  and  Fig.  501,  the  results  for  Cylinder 
3.3d,  Shot  9 appear  more  or  less  reasonable.  While  the  horizontal  force 
indicates  a negative  (upstream)  value  between  about  30  and  100  ms,  it 
is  by  no  means  certain  that  such  is  actually  the  case.  But,  of  course, 
this  is  precisely  the  uncertainty  involved  in  the  present  approach. 

At  later  times,  the  curve  appears  to  oscillate  with  relatively  low  fre- 
quency about  the  dashed  line  indicated  in  Fig.  5-31a-  This  curve  is 
approximately  parallel  to  the  dynamic  pressure  curve  shown  in  Fig.  5 -29, 
and  is  of  the  order  of  magnitude  of  the  dynamic  pressures  acting  at 
this  time.  The  ratio  of  these  two  pressures,  or  "drag  coefficient"^ 
is  about  1.3-  Insofar  as  the  vertical  force  F^  is  concerned,  the 
results  of  Fig.  5.31b  can  only  be  interpreted  as  indicating  a mean 
value  of  zero. 

The  results  for  the  same  cylinder  in  Shot  10,  Fig.  5-32,  are 
similar  to  those  described  above.  The  extent  of  the  negative  portion 
of  Fjj  is  practically  identical  and  the  peak  negative  values  are  in 
the  same  ratio  as  the  incident  overpressure.  Beyond  about  I60  ms,  both 
the  shape  and  magnitude  of  Fj^  seems  to  be  well-defined.  This  portion 
of  the  curve  is  compared  with  the  approximate  dynamic  pressure  curve  in 
Fig.  5 •52a,  and  it  is  seen  that  Fjj  decays  more  rapidly  than  does  the 
dynamic  pressure  curve.  Specifically,  the  "drag  coefficient"  decreases 
from  about  5 ab  I60  ms  to  l.U  at  250  ms.  That  the  drag  pressures  appear 
to  be  unrealistically  high  stands  in  marked  contrast  to  the  fact  that 
initial  peak  forces  are  found  to  be  low.  This  would  tend  to  support 
the  contention  that  initial  peak  forces  cannot  be  obtained  from  the 
strain  data  in  what  might  be  termed  a direct  manner.  However,  an  alter- 
nate interpretation  is  suggested  below.  The  net  vertical  force  again 
indicates  a mean  value  of  zero  with  no  initial  peak  observable. 

Inspection  of  the  original  playback  for  the  cylinder  on  Shot  10 
indicates  a well-defined  frequency  of  about  5 cps,  which  is  also  evident 
in  the  results  for  F^^  and  F^.  The  fact  that  this  frequency  extends 
beyond  the  duration  of  the  blast  suggests  the  possibility  that,  rather 
than  indicating  the  shape  of  the  pseudo-steady-state  loading,  the  results 

^ Quotation  marks  are  used  to  emphasize  the  qualitative  nature  of 
these  results,  both  with  regard  to  the  determination  of  F^  and 
the  particular  choice  of  the  dynamic  pressure  curve. 
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of  Fig.  5.52  simply  represent  a portion  of  this  low-frequency  oscilla- 
tion. In  such  cases,  Fj^  would  then  appear  to  be  oscillating  about 
zero  as  a mean,  and  the  difference  in  "drag  coefficient"  noted  in 
Shots  9 and  10  would  be  due  entirely  to  the  differences  in  amplitude 
of  this  motion.  No  conclusions  can  be  drawn  regarding  these  choices 
of  interpretation. 

With  reference  to  Table  5«5  and  Figs.  5.33a  and  b,  the  results 
for  Cylinder  3»3a,  Shot  10,  deserve  special  mention.  The  strain  data 
for  this  cylinder  were  orginally  discarded,  since  the  strain  calibra- 
tion appeared  to  be  in  error  by  a factor  of  about  10.  However,  the 
results  of  the  graphical  analysis  Indicate,  beyond  about  20  ms,  a hori- 
zontal force  practically  parallel  to  the  dynamic  pressure  curve,  and 
whose  magnitude  yields  an  average  "drag  coefficient"  of  around  0.27; 
this  is  close  to  steady-state  experimental  data  of  about  0.35*  That  the 
peak  force  appears  to  have  been  missed,  of  course,  can  be  attributed 
to  the  complexity  of  the  dynamic  system  noted  earlier.  But  it  is  sur- 
prising that  nearly  all  other  cylinders  showed  a predominant  initial 
peak  — even  if  somewhat  low. 

The  results  obtained  for  the  vertical  force,  Fy,  are  also  of 
interest.  As  can  be  seen  from  Fig.  5 •35b,  Fy  appears  to  Increase 
from  zero  to  a constant  value  at  about  I50  ms;  the  pressure  at  this 
time  is  slightly  in  excess  of  the  dynamic  pressure.  The  trend  indi- 
cated here  for  Fy  appears  to  be  the  only  reasonable  interpretation 
of  the  data  by  the  present  method,  and  stands  in  contrast  to  all  expec- 
tations. It  might  be  noted  that  Cylinder  3*3<i, which  indicated  a zero 
net  vertical  force  in  both  Shots  9 and  10,  was  closer  to  the  ground 
than  Cylinder  3.3a. 

The  data  for  Cylinder  3*3a,  Shot  10,  were  such  that  some 
uncertainty  arose  in  establishing  a contnon  zero  time  in  constructing 
the  horizontal  and  vertical  strain  components.  This  uncertainty 
amounted  to  a relative  time  shift  of  about  5 ms  and,  accordingly,  Fjj 
and  Fy  were  constructed  on  the  basis  of  two  different  zero  positions. 
While  the  resultant  strain  record  differed  to  some  extent  in  phasing, 
the  median  curves  were  consistent  within  the  expected  accuracy  of  the 
method.  The  principal  difference  occurred  in  the  shape  of  Fy  where 
the  initial  rise  time  appeared  to  vary  by  about  20  me,  depending  upon 
the  choice  of  zero  time.  The  results  obtained  for  Fj^  are  Indicated 
in  Fig.  5.53c  where  Fjj  as  determined  from  one  choice  of  zero  time, 
is  superimposed  on  the  strain  record  constructed  according  to  the 
second  choice  in  zero  time.  A comparison  of  the  two  strain  records 
can  be  seen  between  Figs.  5.35a  and  c. 

Cylinder  3.3h  is  the  quarter-scale  model  of  Cylinder  5.5a 
discussed  previously.  The  results  for  Bhot  10  shown  in  Table  5.5  tend 
to  support  the  shape  of  Fy  found  above,  but  are  definitely  not  con- 
sistent with  the  horizontal  force  found  for  Cylinder  5.5a.  As  well  as 
can  be  seen,  Fj^  appears  to  be  oscillating  about  zero  as  a mean,  and 
does  not  indicate  a shape  similar  to  dynamic  pressure  curve.  The  shape 
of  F , however,  is  practically  identical  to  that  found  above,  \dilch 
must  lend  support  to  the  validity  of  this  result. 
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As  discussed  in  Section  ^.2.6  some  vork  was  done  on  the  use 
of  the  ARP  analog  computer  as  a curve-fitting  device  in  connection  with 
the  present  data  reduction  problem.  Preliminary  results  obtained  with 
one  strain  record  on  Cylinder  J.Jb,  Shot  10,  Indicate  a forcing  func- 
tion approximately  parallel  to  the  dynamic  pressure  curves  (see  Fig. 

3.37),  and  nearly  identical  to  the  shape  of  Fjj,  as  found  for  Cylinder 
3.3a^  Shot  10  by  the  present  method.  The  fact  that  this  trend  was  not 
found  by  means  of  graphical  analysis  might  be  due  to  uncertainties  in 
locating  the  common  zero  time,  although,  as  discussed  above,  this  source 
of  error  appeared  to  be  negligible  in  other  cases. 

Cylinder  3 -56  was  the  nearest  to  the  ground  of  the  items  tested. 
As  indicated  in  Table  5*5  and  Fig.  the  shape  of  Fj^  is  consistent 

with  that  found  for  Cylinder  3.5(i  in  both  Shots  9 and  10.  The  initial 
portion  of  Fy  is  not  well-defined  but  beyond  about  80  ms  the  curve 
is  practically  parallel  to  the  dynamic  pressure  curve.  Due  to  the 
qualitative  nature  of  the  analysis,  the  shape  of  F^  found  here  may 
not  be  inconsistent  with  those  found  for  Cylinders  5 *5®  and  b.  Shot  10. 
The  larger  vertical  force  indicated  for  Cylinder  3 -56  is  certainly 
reasonable,  but  shows  a somewhat  smaller  influence  of  ground  proximity 
than  had  been  anticipated. 

5.2.5  Application  of  Transient  Analysis 

A numerical  data- reduct ion  scheme  referred  to  as  transient 
analysis  was  developed  during  the  pretest  planning  phase  of  the  program. 
The  development  and  critical  discussion  of  the  method  is  presented  in 
Appendix  A.  The  present  section  is  concerned  with  the  specific  applica- 
tion of  the  method  to  the  strain  data. 

The  method  of  transient  analysis  was  designed  to  determine  the 
input  or  forcing  function  to  a linear  system  whose  transfer  functions?/ 
are  not  known  but  whose  output  or  response  function  is  known.  In  general, 
for  a multi -degree -of- freedom  system  the  solution  to  this  problem  is  not 
unique;  that  is,  a given  response  in  a particular  mode  of  the  system 
may  result  from  more  than  one  combination  of  forcing  functions  in  the 
various  modes.  In  the  original  development  of  the  method,  it  was 
believed  that  uniqueness  stemmed  from  the  fact  that,  in  the  present 
application,  both  the  unknown  forcing  and  transfer  functions  could  be 
adequately  represented  by  a particular  class  of  analytic  functions. 
Specifically,  the  forcing  function  was  represented  by  a finite  sum  of 
exponential  terms  (i.e.,  "e"  functions),  all  of  whose  exponents  and 
coefficients  were  real,  and  the  transfer  function  was  represented  by 
a similar  series,  all  of  whose  exponents  and  coefficients  occurred  in 
complex  conjugate  pairs.  (This  is  equivalent  to  assuming  a forcing 
function  that  has  no  harmonic  components,  and  a transfer  function  con- 
sisting of  damped  sinusoids.)  Based  on  this  particular  representation, 
the  method  in  brief,  consists  of  curve-fitting  the  known  response  to 


7/  The  transfer  function  refers  to  the  response  of  the  system  in  one 
mode  due  to  a unit  Impulse  applied  in  another  mode  (see  Appendix  A). 
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a similar  exponential  series,  and  then  identifying  those  coefficients 
and  exponents  which  belong  to  the  forcing  function.  (In  the  process 
the  transfer  function  is  determined  also.) 

In  the  present  application  of  the  method,  the  known  response 
was  fitted  to  a five-terra  exponential  series  as  described  in  Appendix  A. 
It  was  hoped  that  the  fifth-order  polynomial  resulting  from  this  curve 
fitting  would  have  three  i-eal  roots  and  two  complex  conjugate  roots. 

If  so,  the  analysis  would  yield  a forcing  function  consisting  of  three 
exponential  terms,  and  a transfer  function  containing  a single  fre- 
quency. However,  as  it  turned  out  in  the  majority  of  cases,  the  analysis 
provided  only  a one-terra  approximation  to  the  forcing  function  and  two 
frequencies.  The  results  obtained  are  shown  in  Table  5*6,  only  a rela- 
tive scale  being  indicated  for  the  magnitude  of  the  forcing  function. 

The  two  cases  for  which  three  terms  in  the  forcing  function 
were  obtained  are  plotted  in  Fig.  5.35;  the  shape  of  the  dynamic  pres- 
sure curve  from  Fig.  5*29  is  indicated  for  comparison.  The  result  for 
Cylinder  5.3c,  Shot  9 (Fig.  5.35a),  shows  no  pronounced  peak,  but  indi- 
cates a rate  of  decay  nearly  identical  to  the  dynamic  pressure  curve. 

(it  will  be  recalled  that  the  results  of  the  previous  section  for  this 
cylinder  indicated  a strong  baseline  shift.)  The  result  for  Cylinder 
5.3e,  Shot  10  (Fig.  5.55h),  is  unrealistic.  In  general,  the  results 
Indicated  in  Table  5.6  are  not  believed  to  be  an  adequate  interpreta- 
tion of  the  strain  data  and,  as  such,  it  must  be  concluded  that  the 
present  application  of  transient  analysis  was  unsatisfactory.  In  an 
attempt  to  understand  Just  why  the  analysis  failed,  a rather  detailed 
post-test  reappraisal  of  the  method  was  undertaken.  Much  of  the  criti- 
cal discussion  presented  in  Appendix  A is  a result  of  this  further 
study. 

The  failure  of  this  first  attempt  can  be  explained  in  terms  of 
the  poor  fit  obtained  for  the  known  response  by  means  of  the  five-terra 
exponential  series.  Whether  it  was  due  to  ccanputational  errors  or 
simply  an  Inadequate  number  of  terms  in  the  approximating  series,  the 
fact  remains  that  the  determination  of  the  forcing  function  was  based 
on  an  analytical  expression  that  was  not  an  adequate  representation  of 
the  measured  strain-time  curve 

As  discussed  in  Appendix  A,  the  valid  application  of  transient 
analysis  is  limited  to  linear  systems  where  the  forcing  functions 
associated  with  each  mode  are  either  zero  or  timewise  proportional. 

When  the  forcing  functions  are  not  all  timewise  proportional,  the 
results  of  the  analysis  probably  cannot  be  given  physical  meaning. 

(This  limitation  of  the  method  was  not  fully  appreciated  at  the  out- 
set.) The  results  obtained  in  the  previous  section  Indicate  that, 

^ The  method  of  transient  analysis  was  applied  to  the  strain  data 
obtained  in  UPSHOT -KNOTHOLE  Project  (Tests  on  the  Loading  of 
Truss  Systems  Common  to  the  Open-Framed  Structures.  WT-723)  with 
equally  unsatisfactory  results.  The  analysis  was  repeated,  based 
on  better  curve-fitting  of  the  data,  but  the  work  was  not  completed 
and  no  conclusive  results  were  obtained.  However,  indications 
were  that  the  analysis  was  capable  of  providing  meaningful  results, 
in  at  least  this  case. 

61 

■seemi  mgfmeTtp 

ONCIASSIFED 


r 


■ ■ 1 

! 

1 

TABLE  5-6  - Results  of  Transient  Analysis;  Coefficients  and  J 

Exponents  of  Forcing  Function,  f(t) 


o,t 

f(t)  * F^e  + F 

Ogt 

2®  + 

a,t 

jC  (t  in  ms) 

Strain 

Numerical  Constants 

Record 

Shot 

m 

wm 

^3 

“2 

^5 

5.5aSl 

10 

1 

1 

-0.017 

5.5aS2 

10 

1 

-0.067 

3.5bSl 

10 

1 

0.4? 

-0.63 

-0.146 

-0.664 

-1.052 

3.3bS2 

10 

1 

-0.013 

3.3cSl 

9 

1 

-0.95 

0.44 

-0.007 

-0.188 

-0.095 

3.3cS2 

9 

1 

1 

-0.147 

3.3dSl 

9 

1 

-0.124 

3.3dS2 

9 

1 

-0.113 

3.3dSl 

1 

0.51 

-0.67 

-0.086 

-0.177 

-0.498 

3.3<iS2 

1 

-0.124 

5.5eSl 

9 

1 

1 

1 

-0.038 

3.5eS2 

9 

1 

-0.088 

3.5eSl 

1 

-0.514 

3.3eS2 

1 

1 

-0.72 

0.90 

-0.080 

-0.505 

1 I 

-0.136 

at  least  initially,  the  Lixe  details  of  the  horizontal  and  vertical 
force  components  differ  considerably.  Hence,  there  is  a strong  indi- 
cation that  at  least  two  modes  of  vibration  are  excited  by  timewise 
dissimilar  forces.  The  fact  that,  initially,  the  method  of  transient 
analysis  may  be  invalid  in  the  present  application,  may  also  explain 
the  unsatisfactory  results  obtained. 


5.2.6  Analog  Computer  Method 


Work  was  initiated  toward  the  end  of  the  post-test  period  on 
use  of  the  ART  Analog  Computer  to  deduce  the  desired  input  functions 
from  the  strain  records.  The  method  of  solution  consists  of  simulating 
a dynamic  system  on  the  computer  which  yields  displacements  similar  to 
the  strain  record  under  consideration,  and  then  adjusting  the  circuit 
parameters  so  that  the  two  displacement-time  curves  are  essentially 
the  same.  The  system  equations  are  taken  to  be  a series  of  linear 
second-order  differential  equations  with  constant  coefficients,  an 
assumption  which  is  also  basic  to  the  method  of  transient  analysis. 

The  forcing  function  may  be  represented  by  a variety  of  analytical  or 
graphical  functions  on  the  computer  but  to-date  only  decaying  exponen- 
tial terms  have  been  uti’ .zed,  again  as  in  the  method  of  transient 
analysis. 

On  first  thought  this  approach  might  seem  rather  hopeless  in 
view  of  the  large  number  of  unknown  quantities  involved;  i.e.,  the 
constants  of  the  system  equations  and  the  parameters  of  the  forcing 
function.  However,  the  method  appears  to  be  quite  practical  and,  while 
not  carried  sufficiently  far  as  to  provide  any  conclusive  results, 
gives  indications  of  being  an  advantageous  approach.  The  feasibility 
of  the  method  in  practice  stems  from  the  fact  that  the  strain  record 
to  be  matched  can  be  plotted  on  graph  paper  and  the  displacement  aris- 
ing from  the  simulated  system  can  be  plotted  automatically  by  means  of 
an  "x-y  recorder''  on  the  same  sheet  of  paper.  This  allows  a visual 
estimate  of  the  accuracy  of  the  matching,  and  permits  much  more  judg- 
ment on  the  part  of  the  operator  in  carrying  out  the  solution  than  is 
possible  using  a numerical  method.  This  advantage,  together  with  the 
brief  time  required  for  varying  parameters  and  obtaining  solutions, 
permits  the  matching  to  be  completed  in  a relatively  short  time. 

A block  diagram  for  a simulated  system  incorporating  two  modes  of 
oscillation  (i.e.,  two  frequencies  in  the  transfer  function)  is  shown 
in  Fig.  5*56.  (in  the  actual  analog  circuit  the  damping  coefficients 
and  the  frequencies  can  be  varied  by  single  adjustments.) 

This  analog  curve-fitting  process  based  on  a two-term  exponen- 
tial representation  of  the  forcing  function  was  applied  to  strain  record 
3.5hS2,  Shot  10.  The  resulting  displacement  curve  and  forcing  function 
are  shown  in  Fig.  5-57;  the  original  strain  record  is  shown  in  Fig. 

5.38  for  comparison.  The  results  shown  in  Fig.  5.37  are  only  of  a 
preliminary  nature  and  do  not  necessarily  represent  the  best  possible 
match  between  the  two  curves  which  can  be  obtained  by  the  present 
method.  Even  so,  a comparison  of  Figs.  5-57  and  5.38  indicates  that 
a reasonable  fit  has  been  obtained.  It  will  be  noticed  that  the  phase 
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of  the  simulated  displacement  differs  somewhat  from  the  actual  record. 
This  can  be  corrected  by  introducing  a third  component  to  the  forcing 
function  having  a finite  value  at  zero  time  (positive  or  negative), 
and  which  decays  to  zero  quite  rapidly.  This  additional  forcing  com- 
ponent would  not  affect  the  displacement  record  after  about  50  “s 
where  the  fit  is  quite  reasonable.  It  would,  of  course,  be  difficult, 
if  not  impractical,  to  match  the  actual  strain  record  identically  over 
the  entire  time  range  of  interest,  since  the  actual  system  may  have 
initially  excited  modes  which  cannot  be  observed  sufficiently  well 
and  thus  cannot  be  simulated. 

The  derived  forcing  function  is  compared  with  the  dynamic  pres- 
sure curve  in  Fig.  5-57,  and  it  is  seen  that  the  wave  form  is  approxi- 
mately the  same  for  both  curves.  It  will  be  recalled  that  the  graphi- 
cal approach  described  in  Section  5*2. U provided  no  conclusive  results 
in  this  case;  if  anything  the  horizontal  force  appeared  to  oscillate 
about  zero  as  a mean. 

The  present  discussion  is  by  no  means  intended  as  a comprehen- 
sive study  of  the  application  of  the  analog  computer  to  the  present 
problem  of  data  reduction,  inasmuch  as  this  technique  was  only  initiated 
toward  the  end  of  the  post-test  period.  However,  this  preliminary  work 
has  indicated  a reasonable  expectation  of  ultimate  success,  and  it  is 
strongly  recommended  that  this  work  be  pursued  further  as  part  of  an 
over-all  study  into  the  data-reduction  problems  associated  with  net 
force  measurement  systems. 

One  advantage  of  the  analog  approach  which  concerns  the  forcing 
function  approximation  should  be  emphasized.  It  is  not  necessary  in 
the  analog  scheme  to  be  restricted  to  an  exponential  approximation  of 
the  forcing  function  (as  is  the  case  in  transient  analysis),  since  it 
is  possible  to  generate  many  other  functions  on  the  canputer.  For 
example,  one  form  of  representation  which  could  be  handled  easily  is 
a series  of  straightline  sections  during  the  initial  loading  period 
followed  by  either  a linear  or  exponential  decay  to  zero.  In  fact, 
existing  load  prediction  schemes  could  be  generated  on  the  computer 
and  the  response  of  the  simulated  system  compared  with  the  actual 
response  data. 
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Normalized  Time,  t = r/U  (Note:  r = 3 f't) 

Fig.  S•^  Normalized  Free-Stream  Id:’essure-Time  Curves 
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CONCLUSIONS  AND  RECCMIENDATIONS 


6.1  CONCLUSIONS 

None  of  the  specific  test  objectives  listed  in  Chapter  2 have 
been  completely  satisfied  as  a result  of  the  test  data  obtained,  and 
in  this  respect  the  program  must  be  termed  unsuccessful.  However, 
some  information  of  value  was  obtained,  and  a great  deal  was  learned 
concerning  the  complexities  of  a field  test  program  such  as  this.  The 
present  section  summarizes  the  results  obtained  and  indicates  the 
extent  to  which  the  various  test  objectives  have  been  satisfied. 

The  test  objective  dealing  with  the  effect  of  shock  strength 
could  not  be  realized  from  the  pressure  data  since,  for  one  reason  or 
another,  the  data  exhibited  a randomness  of  behavior  that  effectively 
obscured  any  such  effect.  In  addition,  little  quantitative  information 
was  obtained  from  the  strain  data. 

Since  the  effect  of  shock  strength  could  not  be  determined,  the 
pressure  data  for  the  first  100  ms  was  averaged  and  further  reduced 
to  yield  what  is  felt  to  be  the  best  possible  estimate  of  the  loading 
on  a cylinder  one  radius  above  the  ground.  While  there  exists  some 
doubt  as  to  the  pressure  scale  ( it  is  believed  to  be  low) , valuable 
information  with  respect  to  the  time  details  of  the  loading  was  obtained. 
In  particular,  the  build  up  and  clearing  times  were  found  to  be  essen- 
tially the  same  as  for  a cylinder  in  free  space,  with  the  exception  of 
the  clearing  time  on  the  front  face  which  was  somewhat  longer.  The 
predicted  clearing  time  on  the  front  face  of  an  isolated  cylinder  is 
between  4 and  5 r/U  time  units;  the  test  results  indicate  a clearing 
time  of  between  8 and  10  r/U  time  units  (r  * radius  of  cylinder, 

U * speed  of  shock  front) . 

The  test  objective  dealing  with  the  effect  of  size  was  to  be 
achieved  by  comparing  the  loading  on  two  geometrically  similar  cylin- 
ders which  differed  in  absolute  size  by  a factor  of  4.  The  pressure 
data  for  these  cylinders  showed  a considerable  and  generally  random 
spread,  and  it  was  not  possible  to  determine  values  of  drag  coefficients 
or  even  a qualitative  trend  in  the  loading  due  to  any  size  effect. 
However,  the  pressure  data  did  indicate  that  the  time  details  of  the 
loading  scale  as  predicted  with  variation  in  size. 
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The  test  objective  dealing  with  the  effect  of  ground  proximity 
could  not  be  realized,  inasmuch  as  four  of  the  six  critical  pressure 
gage  records  (i.e.,  the  gages  on  the  undersurface  of  the  cylinders) 
were  either  lost  or  were  invalid.  In  addition,  the  other  pertinent 
records  are  believed  to  be  poor  with  respect  to  both  magnitude  and 
time  details.  However,  some  Indication  as  to  the  effect  of  ground 
proximity  was  obtained  from  the  strain  data  and  from  the  loading  on 
the  cylinder  one  radius  above  the  ground.  In  the  last-mentioned  load- 
ing a net  upward  vertical  force  of  the  order  of  the  horizontal  force 
was  observed  during  the  drag  phase. 

The  test  objective  dealing  with  the  effect  of  three-dimensionality 
was  to  have  been  achieved  by  comparison  of  the  net  load  as  determined 
by  the  strain  gage  data  with  the  net  (two-dimensional)  load  as  deter- 
mined from  the  pressure  data.  Inasmuch  as  the  strain  data  provided 
virtually  no  quantitative  results,  this  objective  could  not  be  satis- 
fied. However,  some  useful  information  was  obtained  from  the  pressure 
gages  on  the  Project  5 •26.1  vertical  cylindrical  storage  tank,  which 
were  analyzed  in  the  present  report.  These  data  indicated  that  the 
finite  size  of  the  tank  (the  gages  were  close  to  the  roof)  did  not 
affect  the  front  face  loading  appreciably;  the  loading  on  the  back  side 
built  up  somewhat  more  rapidly  than  predicted  for  a two-dimensional 
cylinder,  but  this  is  not  surprising. 

The  strain  gage  data  did  not  provide  anywhere  near  the  infonnation 
desired.  While  there  exists  some  uncertainty  as  to  the  magnitudes  of 
the  recorded  strain,  the  principal  source  of  difficulty  stemmed  from 
an  inability  to  properly  interpret  the  data.  The  cylinder- support 
arrangement  behaved  as  an  extremely  complicated  dynamic  system,  and 
it  was  not  feasible  to  compute  analytically  the  response  of  this  system 
to  a known  input.  There  are  at  least  four  principal  frequencies  evi- 
dent in  most  records,  and  there  is  reason  to  believe  that  motion  of 
the  foundation  in  the  soil  took  place  and  was  picked  up  in  the  strain 
gages.  By  inspection  alone,  the  appearance  of  the  strain  records  was 
such  as  to  yield  only  a very  qualitative  estimate  of  the  forcing 
function. 

The  numerical  data-reduction  scheme,  referred  to  as  transient 
analysis  and  devised  to  interpret  the  strain  data  under  Just  such  cir- 
cumstances, has  proved  unsatisfactory  to  date.  This  is  believed  due 
to  two  main  causes.  First,  a reappraisal  of  the  method  indicates  that 
its  application  to  multi-degree -of -freedom  systems  is  valid  only  where 
the  input  functions  to  the  various  modes  are  either  zero  or  timewise 
proportional.  While  indications  are  that  this  requirement  is  violated 
in  the  present  case,  at  least  during  the  initial  loading  period,  the 
applicability  of  the  method  is  not  sufficiently  well  understood  to 
definitely  rule  out  its  use.  The  failure  of  the  present  application 
of  the  method  probably  can  be  explained  more  readily  in  terms  of  an 
inadequate  representation  of  the  numerical  data  by  means  of  an  exponen- 
tial series  — a step  that  is  crucial  to  the  success  of  the  analysis. 

The  implication  here  is  that  a much  more  refined  Job  of  curve-fitting 
is  necessary  which,  in  effect,  may  increase  the  computational  problems 
enormously.  Therefore,  on  the  basis  to  present  knowledge,  it  is  not 
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possible  to  conclude  definitely’  as  to  the  feasibility  of  the  method 
of  transient  analysis  In  the  present  application.  However,  indications 
now  are  that  it  may  not  be  too  practical  a method  even  though  under 
certain  conditions,  fundamen-tally  sound. 

An  attempt  was  made  to  utilize  the  ARF  analog  ccanputer  in  what 
can  be  termed  a rapid  trlal-and-error  approach  to  the  problem  of  da-ta 
interpretation.  This  work  was  not  completed,  and  while  no  conclusive 
results  were  ob-tained,  the  method  appeared  to  be  a reasonable  one  and 
to  have  good  potentialities. 

It  should  be  emphasized,  however,  that  the  entire  subject  of  da-ta 
reduction  stemming  from  net  force  measurement  systems  deserves  a more 
comprehensive  research  effort  than  could  be  given  to  it  as  a part  of 
the  present  program.  One  cannot  escape  the  fact  that  some  method  of 
data  interpretation  embodying,  at  least  in  principle,  the  philosophy 
of  transient  analysis  is  mandatory,  unless  the  response  of  the  system 
under  consideration  is  such  that  the  desired  information  can  be  ob-tained 
simply  by  inspection. 

There  are,  of  course,  certain  lessons  to  be  learned  from  the 
experiment.  It  now  appears  most  desirable  to  have  obtained  as  much 
information  as  possible  regarding  the  dynamic  systems  being  tested. 

This  would  include  pre-  and  post-test  static  pull  tests  of  the  cylin- 
der (which  would  at  least  have  checked  out  the  strain  gage  system  and 
calibration),  as  well  as  some  sort  of  dynsunic  testing  to  determine  the 
principal  frequencies  and  damping  of  the  system.  Information  of  this 
type  would  have  undoubtedly  aided  in  the  study  and  interpretation  of 
the  strain  data.  For,  if  the  response  of  the  cylinders  to  a known 
input  could  be  obtained  with  sufficient  accuracy,  the  determination  of 
■the  forcing  function  from  the  strain  data  would  have  been  enormously 
simplified. 

In  view  of  the  unsuccessful  application  of  transient  analysis, 
and  for  lack  of  a better  approach,  an  interpre-tation  of  the  strain  da-ta 
was  attempted  on  the  basis  of  inspection  and  simple  graphical  suppres- 
sion of  the  obvious  harmonic  components  of  the  response.  This  approach 
yielded  rather  conflicting  results  which  are  not  without  some  merit  — 
if  only  to  point  out  the  Inadequacy  of  the  method,  or  of  the  strain 
da-ta  itself.  In  one  lone  case  the  horizontal  component  of  applied  load 
substantiated  both  an  approximate  dynamic  pressure  wave  form  and  the 
steady-state  drag  coefficient,  but  there  was  no  indication  of  any 
initial  peak  force.  In  other  cases  the  peak  force  was  evident  and  the 
•wave  form  seemed  reasonable,  but  the  drag  coefficients  were  unrealisti- 
cally high.  In  still  other  cases,  the  rate  of  decay  of  the  force 
differed  significantly  from  the  approximate  dynamic  pressure  curve. 

The  data  from  two  cylinders  indicated  a vertical  force  that  rose  gradu- 
ally to  a value  in  excess  of  the  dynamic  pressure  at  that  time;  cylin- 
ders closer  to  the  ground  indicated  a vertical  force  that  oscillated 
about  zero  as  a mean. 
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6.2  RECOMMENDATIONS 


While  the  present  test  has  been  largely  unsuccessful,  it  is 
believed  that  the  original  objectives  are  still  valid  ones.  In  view 
of  the  support  lent  the  scaling  of  time  details  and,  in  some  cases  , 

pressures,  it  is  recommended  that  additional  experimental  investiga- 
tions of  the  blast  loading  on  cylinders  be  conducted  on  a model  scale 
in  a shock  tube  utilizing  pressure,  interferometric,  and  direct  force 
measurements. 

Much  additional  work  remains  to  be  done  with  respect  to  net 
force  measurement  systems,  both  as  to  physical  design  and  interpreta- 
tion of  data.  Despite  the  present  unsuccessful  utilization  of  the  l 

cylinder  as  a dynamometer  with  strain  output,  this  approach  has  obvious 
advantages.^  It  is  recommended,  therefore,  that  a comprehensive 
study  be  initiated  on  the  data-reduction  problems  associated  with  net 
force  measurement  systems.  The  objectives  of  this  study  would  be  to 

(a)  devise  methods  of  analysis,  (b)  specify  the  class  of  physical  ; 

systems  whose  response  will  successfully  permit  the  Inverse  solution,  i 

and  (c)  determine  the  type  and  quality  of  the  necessary  data  to  permit 
this  solution.  It  is  felt  that  particular  consideration  should  be 

given  to  the  analog  computer  technique  of  data  analysis  discussed  in  ■; 

this  report . I 

i 

ii 


"y  Of  course,  other  schemes  can  be  thought  of  for  the  present  test 
items  (e.g.,  utilization  of  the  cylinder  as  a ballistic  pendulum 
with  accelerometer  output) . However,  the  present  dynamometer 
arrangement  would  be  applicable  to  a wide  variety  of  items. 


94 


APPENDIX  A 


TRANSIENT  ANALYSIS 


A.l  INTRODUCTION 


In  the  present  application,  the  prime  data  (i.e.,  the  measured 
strains)  do  not  represent  in  a direct  sense  the  desired  information  (i.e. 
the  net  forces).  In  fact,  the  data  obtained  represent  the  response  of 
a rather  complicated  dynamic  system  to  an  input  which  is  itself  the 
desired  information.  The  dynamic  system  consists  of  aerodynsunic, 
elastic,  plastic,  electronic,  thermal  and  possibly  stochastic  components 
coupled  in  many  and  varied  fashions.  The  problem  at  hand  is  to  remove 
from  the  prime  data  the  influence  or  contributions  of  these  various 
components  and  of  the  system  itself,  leaving  the  desired  information. 
Since  it  is  difficult  (probably  impossible)  to  assess  the  exact  contri- 
bution or  influence  that  the  dynamic  system  contributes  to  the  prime 
data,  a rational  approach  to  the  problem  is  to  attempt  to  assess  the 
contribution  or  influence  of  those  components  which  could  be  expected 
to  have  the  most  pronounced  effects.  Due  to  lack  of  sufficient  know- 
ledge of  the  test  conditions  and  the  reliability  and  characteristics 
of  the  various  recording  and  data -reduct ion  Instruments,  the  magnitude 
of  the  contributions  of  some  of  the  above  components  cannot  even  be 
estimated.  Nonetheless,  it  has  been  assumed  that  the  major  influence 
on  the  prime  data  stems  from  an  elastic  response  of  the  cylinder-support 
system . 

The  following  analysis,  referred  to  as  "transient  analysis,"  is 
designed  to  remove  the  influence  of  the  elastic  components  in  the  prime 
data,  which  is  undoubtedly  strongest,  regardless  of  other  effects. 

The  analysis  as  applied  here,  represents  the  function  or  input  required 
to  produce  the  prime  data  as  the  response  of  an  elastic  model. 

A. 2 DERIVATION  OF  THE  BASIC  EQUATION 


For  the  present  application,  the  elastic  model  is  assumed  to  be 
linear  and  possess  several  (but  a finite  number  of)  degrees  of  freedom. 
Its  motion  is  then  determined  by  a system  of  linear  differential 
equations  with  constant  coefficients.  One  such  equation  for  each 
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degree  of  freedom  occurs.  Some  of  the  possible  degrees  of  freedom 
which  might  be  effective  are  the  rigid-body  translation  and  rotation 
of  the  foundation,  rigid-body  motions  of  the  cylinder,  bend  ng  and 
axiaj.  vibrations  of  the  cylinder  supports,  bending  of  the  cylinder, 
torsional  deflections  of  the  cylinder,  or  vibrations  of  the  cylinder 
as  a shell.  Regardless  of  the  number  of  degrees  of  freedom  present, 
the  equations  of  elastic  motion  may  be  written  in  the  form: 


j - 1,  2,  ...  N 


where  the  b^^j,  and  mj  are  real  constants,  and  fj(t)  is  the  forc- 

ing term  for  the  Jth  degree  of  freedom.  Xj(t)  is  the  response  of  the 
Jth  degree  of  freedom  represented  by  the  Jth  mass,  mj.  Equation  A.l 
is  the  most  general  system  of  equations  for  N degrees  of  freedom  with- 
out inertial  coupling.  The  constants  a^^j  and  b^j  are  the  damping 
and  elastic  coupling  coefficients,  respectively,  some  of  which  may  be 
zero. 


For  the  present  application,  the  forcing  terms  fj(t)  arise 
primarily  from  the  pressure  acting  on  the  cylinder.  It  is  reasonable 
to  assume  that  motion  of  the  foundation,  if  it  occurs,  can  result  only 
from  coupling  with  the  cylinder  motion.  Thus,  only  those  of  Equation 
A.l  corresponding  to  masses  ra.  associated  with  motion  of  the  cylinder 
are  nonhomogeneous  (i.e.,  fj(tj  4 O)- 

The  masses  mj  associated  with  the  motion  of  the  cylinder  are 
of  two  types;  namely,  total  mass  corresponding  to  rjgid-body  motions, 
or  generalized  modal  masses  for  the  bending  or  shell  vibrations.  The 
forcing  terms  for  each  of  these  result  from  either  the  integrated 
pressure  over  the  cylinder  alone  or  the  integrated  pressure  times 
a modal  deflection  curve  over  the  cylinder.  If  the  pressure  is  assumed 
to  be  nearly  constant  in  the  axial  direction  along  the  cylinder,  these 
forcing  terms  will  be  essentially  proportional  to  one  another  at  each 
instant  of  time  so  that  their  time  histories  will  be  similar.  In  effect 
Equation  A.l  will  contain  the  same  rlghthand  members  up  to  a constant 
factor,  which  may  be  zero. 

In  this  case.  Equation  A.l  becomes. 


(t)  + 


J = 1,  2 ...  N 


(A. 2) 
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It  is  well-known  (see  Vibration  Problems  In  Engineering  by  S.  Timoshenko) 


that  the  response  Xj(t),  Equation  A.^ subject  to  homogeneous  initial 
conditions,  can  be  given  by  a sum  of  Duhamel-type  integrals  as: 


Xj(t) 


M t 

i I 


(t  - T)i  T 


(A. 5) 


where  the  hj^j  are  referred  to  as  transfer  functions  and,  for  a fixed 
represent  the  impulsive  response  in  the  Jth  degree  of  freedom 
to  a unit  impulse  applied  in  the  kth  degree.  In  other  words,  hjjj(t) 
is  the  solution,  Xj(t),  of  the  corresponding  homogeneous  system  of 
equations 


mjXj(t)  + ° 

J « 1,  2,  ...  N 

subjected  to  the  inhomogeneous  initial  conditions 
Xj(0)  - Xj(0)  - 0 J ^ k 

Xj^(O)  - 0,  Xj^(O)  . 1 

Such  a solution  Xj(t)  ■ ^kjt'*')  known  to  be  a linear  combination  of 
exponential  functions  (in  general,  complex)^  of  the  form  (see  Vibra- 
tion Problems  in  Engineering  by  S.  Timoshenko). 


^ k1 


(A. 5) 


where  Xjj,(n  ■ 0,  tl,  ...  ±N)  are  the  2N  roots  of  the  characteristic 
equation  (system  of  equations) 


^ i«l  ''  i-1  ^ 


J « 1,  2,  ...  N 

If  Equation  A.^  is  substituted  in  Equation  A.J,  there  results: 
N N ^ (t  - t) 

V*’- 1 


Ca.6) 


It  is  assumed  here  that  the  2N  characteristic  roots  are  distinct. 


(A. 6a) 


If,  in  addition,  the  fi^(t)  are  timewiae  proportional  so  that 

fj^(t)  - Cj^f(t) 


as  in  Equation  A. 2,  then.  Equation  A. 6 n»y  be  written, 


Xj(t)  . 

Xj(t)  . 
Xj(t)  . 


J 


N N h.1  Jt  - r) 

e ° d T 

n 


f I 2 

J fc.1  n«-N 

o 

J f(T)  I e 

V n«-N 


^(t)  hj  (t  - t)  dr 


d T 


(A.7) 


(A.  7a) 


It  is  seen  that  the  quantity  hjCt-  t)  in  Equation  A.7  is  precisely 
the  form  of  the  transfer  functions  hjjj(t)  in  Equation  A. 5 (independent 
of  k,  however)  and  is  a linear  combination  of  exponential  functions 
involving  only  the  characteristic  roots  of  the  system.  Equation  A.7 
is  the  desired  basic  equation  and  may  be  used  to  obtain  the  unknown 
forcing  function  f(t)  from  a known  response  Xj(t).  The  quantity  f(t) 
actually  represents  only  the  shape  of  the  forcing  function,  since  It  is 
defined  up  to  the  constant  factors  C*  of  Equation  A. 2. 

The  problem  of  uniqueness  and  the  effect  of  "errors"  in  the 
response  curve  Xj(t)  as  they  affect  the  solution  for  f (t)  obtained, 
will  be  discussed  at  the  appropriate  steps  in  the  subsequent  analysis. 

It  will  be  seen  that  the  uniqueness  problem  essentic^ly  reduces  to 
a problem  in  curve-fitting  and  may  be  treated  in  several  ways. 

A. 3 SOLUTION  OF  THE  BASIC  EQUATION 


Equation  A.7  can  be  solved  by  a variety  of  methods,  none  of 
which,  however,  leads  to  a unique  solution  for  the  forcing  function 
f(t).  For  the  present  application  it  is  assumed  that  f(t)  can  be 
represented  by  exponential  functions  in  the  form: 
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f(t) 


F^e  + 


V 


(A.8) 


1>1 


P^e 


ait 


where  the  F^  and  a.  are  constants.  The  number  (K)  of  terns  in  Equation 
A.8  will  be^detennlned  by  the  complexity  of  the  actual  forcing  function 
or,  at  least,  the  complexity  required  of  f(t)  in  order  to  produce  the 
known  response  X(t).  It  should  be  mentioned  that  the  fonn  of  Equation 
a.8  in  no  way  restricts  the  generality  of  application. 

If  Equation  A.8  is  substituted  into  Equation  A. 7,  and  Equation 
A. 7a  is  used,  there  results  upon  integration: 


X^(t) 


K 

i-1 


Aie 


V 


N 

y B e 
n»-N 


(4-9) 


where. 


N 

^ n 

t^-N  i n 


"1 

i»l  ®i  ■ ^n 


(A. 10) 


It  is  assumed  here  that  ^ "x  since  other- 

wise a resonance  condition  would  result,  which  certainly  does  not 
happen  in  the  physical  case  considered  here.  Equation  A. 9 gives  Xj(t) 
as  a sum  of  exponential  functions,  and  the  present  problem  reduces 
to  that  of  fitting  such  a sum  to  the  prime  data  representing  Xj(t) 

for  a discrete  set  of  tiroes,  t. 

A curve-fitting  scheme  results  in  the  constants,  Ai,  a^,  B^^, 
and  Xn  (real  and  complex,  in  general)  which,  in  turn,  may  be  used  to 
compute  the  desired  F^  through  Equation  A. 10.  Actually,  Equation 
A. 10  will  also  yield  the  which,  however,  are  not  of  direct  concern 
in  the  present  application.  With  the  a^,  B„,  and  know, 

Equation  A. 10  constitutes  a system  of  2N+K  equations  in  the  2N+K 
unknowns,  and  F^.  With  the  a ^ and  F^  determined,  the  desired 

forcing  function  is  then  given  by  Equation  A.8.  Actually,  a little 
more  information  is  assumed  here;  namely,  that  not  only  the  constants 
a.  Ai,  Bn,  and  Xj,  are  obtained  from  the  curve-fitting,  but  that 
the  sets  of  A^'s  and  a .'s  can  be  distinguished  from  the  set  of 
B *8  and  X 's.  Otherwise,  the  structure  of  Equation  A. 10  would  not 
be  determined  so  that,  a unique  forcing  function  could  not  be  constructed 
through  Equation  A.8.  In  the  present  application,  the  distinction 
between  the  two  classes  of  tenns  in  Equation  A. 9 can  be  made  through 
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the  exponents 


themselves.  Since  Equation  A.l  can  be 


assumed  to  be  stable,  none  of  the  characteristic  roots  Xq  can  be 
pure  real.  On  the  other  hand,  It  Is  assumed  that  the  forcing  function 
f(t)  contains  no  harmonic  components  so  that  a realistic  curve-fitting 
scheme  should  yield  only  real  exponentials  Oj  for  f(t)  In  Equation 
A. 8.  In  brief,  the  nonreal  exponents  obtained  from  curve-fitting  are 
associated  with  the  transfer  function 

b,(t)  . ■£  iSj  e ° 

and  give  the  \dille  any  real  exponents  obtained  are  associated 
with  the  forcing  function 

K a.t 

f'(t)  - i;  P.e  ^ . 

1-1 

The  appropriateness  of  this  choice  of  association  can  be  asserted  only 
after  sane  study  Is  made  of  the  effect  of  errors,  etc.,  of  the  prime 
data  on  the  computed  exponentials.  It  would  be  desirable,  of  course, 
to  obtain  the  natural  frequencies  of  the  system  either  analytically 
or  experimentally  Independent  of  the  above  analysis,  since  then  they 
could  be  Identified  in  Equation  A. 9 and  f(t)  could  then  contain 
harmonics  without  affecting  the  above  solution. 

A. 4 CURVE-FITTING 

Equation  A. 9 may  be  written  in  the  form: 


M y^t 
X(t)  - 2 C e 
k»l  * 


(A. 11) 


where  the  Cj^  and  are  in  general  complex  numbers,  and  where  the 
subscript  j has  been  dropped,  since  the  only  concern  is  with  the 

response  of  but  one  component  of  the  system.  Curve-fitting  consists 


of  finding  those  values  for 


such  that 


as  given  by 


Equation  A. 11  best  fits  the  given  prime  data. 

Method  A.  The  most  direct  method  of  curve-fitting  is  accomplished 
by  visual  inspection  of  plots  of  the  prime  data  itself.  A discussion 
of  this  method  and  the  results  of  its  use  in  the  present  application 
are  presented  in  Section  5.1.2. 

Method  B.  A second  method  consists  of  rather  standard  technique 
from  the  theory  contained  In  the  calculus  of  finite  differences  (see 
The  Calculus  of  Observations  by  Whittaker,  E.  and  Robinson,  A.).  It  is 
assumed  now  that  the  prime  data  are  given  in  equal  inteirval  tabular 
form. 
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Let  the  interval  of  tabulation  of  the  function  X(t)  be  A , and 
designate  X(t)  at  the  time  when 


t = n A 


n = 0,  1,  . . . , P 


as  Xjj.  The  assumed  form  of  X(t)  (Equation  A. 11)  may  then  be  written 


(A. 12) 


where 


y.  A 


X,  = e 


(A. 12a) 


This  method  of  curve-fitting  rests  upon  the  observation  that  Equation 
A. 12  is  the  solution  to  a linear  finite  difference  equation  with 
constant  coefficients.  This  finite  difference  equation  may  be  written 
for  each  consecutive  group  of  (L  + l)  tabulated  values  of  X^  as 


1?0  ° 


(A. 15) 


where  = 1 and  the  remaining  D's  are  to  be  determined.  If  (P  + 1) 
is  the  total  number  of  tabulated  values  of  Xjj,  Equation  A. 15  may  be 
written  (P  - L + l)  times,  each  equation  corresponding  to  each  of  the 
values  of  n in  the  range 

n s 0,  1,  2,  •••,  (p  - L) 

Consider  the  resulting  set  of  simultaneous  equations 


^0^0  * ^1^1  ^ • • • + ®L-1  ^L-1  “ ” 


DqX^  + ' 


^O^P-L  ^ ^l^P-L  ^ ^ 


(A.llf) 
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where  the  X ’s  are  known  from  the  tabulated  function  X(t).  This  set 
of  simultaneous  equations  consists  of  (P  - L + l)  equations  in  the 
(L  + 1)  unknown  D’s.  A solution  to  the  system  of  equations  is  possi- 
ble if  P^2L.  When  L,  the  number  of  exponential  terms  used  to 
approximate  X(t),  is  chosen  equal  to  P/2,  the  system  has  a unique 
solution  if  the  determinant  of  the  coefficients  of  Xn  1® 

Normally,  this  condition  exists  and  the  unknown  D's  are  determined  by 
ordinary  methods  (see  The  Calculus  of  Observations  by  Whittaker, E. 
and  Robinson,  A.)  and  are  used  to  form  the  characteristic  equation 
which  is  the  polynomial 

0-D  +D  X + dX^+  ,.,+  D-.X^^  + X^  (a.  15) 

0 1 2 1>-1 

The  roots  of  Equation  A. 15  are  the  X^'s  referred  to  in  Equation  A. 12 
and,  hence,  may  be  used  to  determine  the  Xk'®  turn,  are 

used  to  form  the  approximate  representation  of  X(t). 

The  actual  determination  of  the  Xy^'s  is  considered  now,  since 
the  roots  of  Equation  A. 15  may  be  positive,  negative,  or  complex. 
Regardless  of  the  type  of  root,  X ^ can  be  written  as 


(i.^) 

where  8 is  equal  to  zero  or  pi  ( TT ) if  Xj^  is  positive  or  negative, 
respectively.  If  Xy^  is  a complex  root,  0 is  the  Cauchy  principle 
angle  and  I'xes  between  - v and  TT  . I^  Xj^  is  positive. 


(A. 16a) 


If  X^  is  negative. 


X 


k 


log  (-x^)  ♦ 


1 TT 

A 


(A. 16b) 


In  the  latter  case  the  imaginary  part  of  Xj^  is  dropp^,  since  it 
introduces  an  extraneous  frequency  whose  period  is  2Za.  The  con- 
clusion that  this  frequency  is  extraneous  is  based  upon  the  observa 
tion  that  it  is  dependent  only  upon  the  choice  of  the  tabulation 


I 
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interval  and,  therefore,  may  be  changed  to  any  value  by  a suitable 
choice  of  interval. 2/  Thus,  for  all  real  values  of 


(A. 17a) 


If  Xjj  Is  complex,  it  must  be  one  of  a conjugate  pair  and  both  pairs 
are  contained  in 


(A. 17b) 


Having  obtained  the  y^'s  from  Equations  A. 17a  and  A. 17b,  the 
Cjt's  may  be  obtained  from  Equation  A. 12  in  combination  vith  the  initial 
conditions.  These  conditions  are  satisfied  if 


M 

VmI 


and 


M 


1^1 


Vk 


(A.  18) 


(A.  19) 


2/  An  alternate  method  of  dropping  the  frequency  induced  by  the  nega- 
tive roots  is  to  consider  only  the  even  terms  in  Equation  A. 12. 
Equation  A. 12  becomes 

k«l 

And  Equation  A. 12a  becomes 


from  Tdilch 


and  X(t)  again  has  the  form  of  Equation  A. 12  where  the  extraneous 
frequency  corresponding  to  IT /A  present. 
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Since  Equation  A. 12  is  to  be  satisfied  for  each  tabulated  value  of 
X(t),  there  results  the  following  P equations  in  (L  - 2)  unknowns, 


C ^ X ^ + Cg  ^ 2 ^ • • • + Cj^  ^ L * ^1 
X ^ + Cg  Xg  + . . . + ^ L * ^2 


*^1^1  ^2^2"^ 


(a. 20) 


where  Equations  A.l8  and  A. 19  are  to  be  used  to  eliminate  two  of  the 
Cjj's.  The  above  system  of  equations  is  overdeterminate  and  usually 
inconsistent,  since  a finite  number  of  terms  is  used  to  approximate 
X(t).  The  method  of  least  squares  is  used  to  detennine  the  most 
probable  set  of  (L  - 2)  equations  in  the  (L  - 2)  unknowns  (see  Numerical 
Mathematical  Analysis,  by  J.  B.  Scarborough)  and  then  this  set  of 


equations  is  solved  for  the 


Method  C.  A third  method  of  curve-fitting  utilizes  an  andlog 
computer,  and  represents  a rapid  trial-and-error  process.  A discussion 
of  this  method  is  presented  in  Section  5-l»5« 

A. 5 DISCUSSION  OF  APPLICATION  OF  TRANSTEW?  ANALYSIS  (METHOD  B) 


The  application  of  transient  analysis  as  described  in  Method  B 
consists  of  the  following  steps: 

a)  The  frequencies  and  exponents  occurring  in  the  expansion 
of  the  response  curves  (prime  data)  are  determined  by 
analogy  with  finite  difference  equations.  The  prime  data 
is  interpreted  as  representing  the  solution  of  a system 
of  linear  finite  difference  equations.  The  character- 
istic equation  (polynomial)  for  this  system  of  difference 
equations  is  constructed  using  the  prime  data  and  solved 
to  obtain  the  characteristic  roots.  These  roots,  in  turn, 
are  taken  as  the  desired  frequencies  and  exponents  of 

the  response  curve  expansion. 

b)  Having  obtained  the  frequencies  and  exponents  of  the 
response  curve  expansion,  the  suoplltudes  associated  with 
each  of  the  exponential  terms  are  then  determined  by 

a straightforward  least  squares  fit  to  the  prime  data. 

c)  With  the  expansion  of  the  response  curve  determined,  the 
forcing  function  is  then  associated  with  the  reeil  expo- 
nential terms  and  is  obtained  through  Equation  A. 10. 

2/  If  L<P/2,  the  method  of  least  squares  is  used  to  reduce  the 
equations  to  a determinate  number. 
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I The  success  of  the  application  of  transient  analysis  to  the 

I predicting  of  actual  forcing  functions  hinges  upon  several  reasonable 

' but  crucial  assmprtlons.  First,  and  of  utoost  importance,  is  the 

" ' assumption  that  the  dynamic  system  poroducing  the  response  curves  be  f 

essentially  linear.  Just  \diat  magnitude  of  nonlinear  effects  such  as  j 

stochastic  errors,  or  baseline  shifts  are  tolerable  is  not  Imoun  at 

• present,  indeed  represents  a major  area  of  required  investigation  | 

before  any  net  force  measurement  technique  of  the  present  type  can  be  ± 

evaluated.  The  effect  of  errors  (known  or  unknown)  in  the  prime  data 

affect  the  analysis  most  critically  in  step  a)  above.  In  fact,  if  the 

frequencies  and  exponents  determined  in  step  a)  are  realistic,  the  i 

success  of  the  entire  analysis  appears  to  rest  only  upon  the  Identifica:- 

tlon  made  in  step  c),  idvLch  will  be  discussed  presently.  One  criticism 

of  the  procedtire  in  step  a)  for  determining  the  all-important  fre-  | 

quencies  and  exponents  appears  immediately  as  the  result  of  using  the  | 

prime  data  Itself.  Since  these  data  in  the  present  application  are  | 

loiown  to  contain  various  types  of  errors,  it  would  be  desirable  to  I 

incorporate  some  datar-smoothing  process  prior  to  the  application  of 

transient  analysis.  There  are  several  types  of  datar-smoothing  techniques 

available  with  various  degrees  of  reliability  and  meaningfulnoss  to  the 

present  application,  but  the  most  promising  appears  to  be  the  use  of  the  y 

autocorrelation  function.  For  the  latter  application,  step  a)  would  use  I 

• the  autocorrelation  function  defined  as  I 


T 

X(t)  = P X(T)X(t  + t)  d T (A. 21) 

0 0 

instead  of  the  prime  data,  X(t),  itself.^ 

The  autocorrelation  function  ir(t)  contains  the  same  frequencies 
and  exponentials  as  does  the  prime  data,  X(t),  and  smooths  the  data 
by  eliminating  stochastic  and/or  bias  errors.  This  smoothing  technique 
is  being  tried  and  tested  at  present  in  other  and  similar  applications, 
(see  Operation  TEAPOT  Project  3.2,  Study  of  Drag  Loading  on  Structures 
in  and  out  of  the  Precursor  Zone,  WT-1124yT  Of  course,  the  analogy 
technique  applied  in  step  a)  does  itself  perform  some  sort  of  averaging 
on  the  data  fed  in.  Just  what  the  nature  of  this  averaging  is,  is  not 
known  at  present.  A perturbation  study  of  the  effect  on  the  exponents 
and  frequencies  as  determined  by  step  a)  is  necessary  before  a realis- 
tic evaluation  of  transient  analysis  can  be  made.  Because  of  the 
complications  involved,  such  a study  would,  no  doubt,  have  to  be  car- 
ried on  in  the  sphere  of  numerical  experimentation. 

Step  b)  determines  a "best  fit"  of  the  prime  data.  The  fit 
obtained  itself  is,  in  many  respects,  a measure  of  the  reliability  of 
transient  analysis  in  the  present  application.  However,  a realistic 
criteria  for  "goodness"  of  fit  has  not  yet  been  established;  at  least, 
not  in  terms  of  the  reliability  or  meaningfulness  of  the  forcing  function 


y The  fixed  interval  T is  taken  approximately  as  representing  the 
extent  of  the  data  available. 


80  obtained.  A point-by-point  check  on  the  fit  may  not  be  too  meaning- 
ful, since  the  analysis,  to  be  successful,  must,  of  necessity,  ignore 
"errors"  in  the  prime  data  and  pick  out  only  pertinent  information. 
However,  the  method  of  least  squares  as  applied  in  step  b)  appears, 
at  present,  to  have  the  most  reliable  characteristics  for  curve-fitting 
as  applied  in  this  application. 

The  procedure  of  step  c)  is  based  on  the  assumption  that  the 
forcing  function  is  represented,  or  "sufficiently  well"  represented, 
by  a sum  of  real  exponentials.  Just  what  is  "sufficiently  well"  to 
assure  success  of  transient  analysis  has  not  yet  been  determined.  The 
limits  of  reliability  could  most  likely  (and  probably  have  to)  be  deter- 
mined by  numerical  experimentation,  also. 

Finally,  the  assumption  that  all  forcing  functions  are  timewise 
proportional  needs  to  be  considered.  The  first  question  to  be  answered 
is  whether  the  forcing  functions  considered  are  essentially  propor- 
tional or  not.  The  second  question  to  be  answered  is  what  is  "suffi- 
ciently proportional"  in  terms  of  success  of  the  application  of  tran- 
sient analysis.  Thirdly  (since  the  first  two  questions  have  not  as  yet 
been  answered),  what  happens  to  transient  analysis  if  the  forcing 
functions  are  not  essentially  timewise  proportional?  In  jjarticular, 
does  the  forcing  function  obtained  by  applying  transient  analysis  in 
this  latter  case  have  any  significance  or  yield  any  real  information? 

Of  these  three  questions,  only  the  third  can  be  answered  (to  some 
degree  at  least)  at  the  present  time. 

In  answer  to  the  question  of  what  happens  if  the  forcing  functions 
are  not  timewise  proportional,  the  analysis  of  the  preceding  sections 
is  repeated  in  brief  as  follows: 

Equation  A. 6 is  written  in  the  form 


X(t) 


N t N X„(t-r) 

2 f 2 

k*l  J-  n.-N 


(A. 22) 


where  the  index  j has  been  dropped,  since  the  only  concern  is  the 
response  in  a single  degree  of  freedom.  The  forcing  function  fk(t) 
is  written  in  the  form 


K . a.t 

flj(t)  - 2 (A. 23) 

Isl 

where  the  a range  over  all  exponentials  appearing  in  any  of  the 
degrees  of  freedom,  so  that  some  of  the  may  be  zero.  Vlith 

Equation  A. 23  substituted  in  Equation  A. 22  and  the  Integration  per- 
formed, there  results  the  equation: 
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(A.2U) 


which  is  precisely  the  same  form  as  Equation  A. 9.  However,  the 
and  Bq  are  related  to  the  forcing  coefficients  and  the  transfer 

function  coefficients  0^^  by 


N N 


_k 
^i  °n 


2 2 5-Vr 

k^l  n— N i ''n 


^ “n 


N K 

B.  = 2 I 

fa=l  i*l 

which  are  the  countei-parts  of  Equations  A. 10.  The  relationships 
between  Equations  A. 10  and  A. 25  result  in  (dropping  the  index  j) 


(A. 25) 
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T F 

W " " ^ni 


(A. 26) 


n - tl,  ^2,  ...  *N 

Isl,  2,  ...K 

Hence,  the  curve-fitting  of  step  b)  may  be  applied  to  Equation  A. 24 
resulting  in  the  A^,  a^,  Bq  and  Xq  which,  in  turn,  are  used  in 
Equation  A. 10  to  determine  the  and  F^.  However,  step  c)  Is  not, 
in  general,  valid  since  the  F^  are  not  the  actual  forcing  function 
coefficients.  In  fact,  there  is  no  longer  a single  forcing  function 
in  question,  but  one  for  each  degree  of  freedom.  The  actual  forcing 
function  coefficients  F^  are  related  to  the  and  F^  by 

Equation  A. 26.  However,  in  Equation  A. 26  only  the  /9n  and  F^ 
are  known,  while  both  sets  of  Ff  and  are  unknown.  Since  there 

are  (2N)  x K equations  in  Equation  A. 26  and  (2N)  x K x N unknowns. 
Equation  A. 26  does  not,  in  general,  determine  the  actual  forcing 
function  coefficients  F^  uniquely  and  an  infinitude  of  solution  may 
be  constructed.  It  should  be  mentioned,  however,  that  some  informa- 
tion about  the  forcing  functions  has  already  been  determined  in 
step  b) ; namely,  aU  the  exponents  a ^ appearing  in  ^ of  forcing 
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functions  have  been  determined,  and  Equation  A. 26  represents 

a number  (even  though,  in  general,  not  enough),  of  necessary  conditions 
on  the  coefficients  . Certain  auxiliary  information  about  either 
or  both  the  dynamic  system  or  forcing  functions  obtained  analytically 
or  experimentally  together  with  Equation  A. 26  can,  in  certain  instances, 
determine  a unique  set  of  coefficients  f]^  . For  example,  if  curve- 
fitting results  in  but  one  real  exponential,  transient  analysis  yields 
a unique  forcing  function  (in  a trial  way,  of  course). 

It  is  interesting  to  see  how  the  assumption  of  timewise  propor- 
tionality of  forcing  functions,  together  with  Equation  A. 2^  yields  a 
unique  solution  (as  it  must,  by  previous  considerations). 

The  assumption  of  proportionality  is  equivalent  to  Equations 
k.2'k,  A. 6a  and  A. 8: 

F^  a each  i and  k (A. 27) 

so  that  Equation  A. 26  may  be  written 


fel 


a ^ 

kin 


^ k 
k»l 


^1  “ 

n = ti,  ±2,  ...  ±N 
i = 1,  2,  ...  K 


(A. 28) 


which  identifies  the  F^^,  as  well  as  the 
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n 
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k 

n 


of  Equation  A.7a^  uniquely,  since  the  right  members  of  Equation  A. 28 
are  known.  In  effect  the  assumption  of  proportionality  uncouples 
Equation  A. 26  (which  is  generally  coupled  in  the  independent  indices 
i and  n)  into  two  independent  sets  of  equations 


and 


F^  a F^  ; i = 1,  2,  . . . , K 


y C,  ; n = tl,  ±2,  ...  ±N 

k=l 


^ Dropping  the  index  J. 
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Less  draBtic  assumptions  regarding  the  forcing  functions  or 
transfer  functions  could  be  expected  to  yield  additional  information 
about  the  forcing  functions.  For  example,  if  Equatio’-  £ .26  is  summed 
with  respect  to  i,  there  results 
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k*l 
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_i=l 

(A.29) 


The  inner  sum 


F 


k 
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is  clearly  (see  Equation  A. 23)  the  quantity  fjj(O)  which  is  the  magni 
tude  of  the  forcing  function  in  the  kth  component  at  time  t ® 0. 
Hence,  Equation  A.29  “ay  be  written: 
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E \i0)  a 

k=:l 
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(A.30) 


n « tl,  ±2,  ....  ±N 

It  is  possible  that  at  least  the  relative  magnitudes  (if  not 
the  absolute  magnitudes)  of  the  fk(0)  can  be  estimated  from  theory, 
etc.,  so  that  Equation  A.30  could  be  used  to  yield  additional  inde- 
pendent conditions  on  the  unknowns  of  Equation  A. 26.  If  the  absolute 
magnitudes  of  the  f^(0)  are  knoira.  Equation  A.30  yields  2N  additional 
equations.  If  only  their  relative  magnitudes  are  known,  2N-1  addi- 
tional equations  result.  , 

Another  possibility  is  that  some  of  the  F*  may  be  known  from 
other  considerations,  in  effect  reducing  the  number  of  unknowns  in 
Equation  A. 26.  For  exaunple,  some  degrees  of  freedom  may  not  admit 
a forcing  function  and  are  known  to  be  excited  only  through  coupling 
with  other  degrees.  Hence,  the  corresponding  F^  are  zero.  Certain 
timewise  proportionality  conditions  may  also  exist  between  certain 
degrees  of  freedom,  such  as  the  set  of  all  bending  modes,  further 
reducing  the  number  of  unknowns. 
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10  The  Quartermaster  General,  D/A,  Washington  2^, 

D.C.  ATTN:  Besearch  and  DevelopMnt  Dlv. 

11*  15  Chief  of  Engineers,  D/a,  Washington  25,  D.C.  ATW; 
ENGRB 

lo  Chief  of  Transportation,  *llltary  Planning  and  Intel- 
ligence Dlv.,  Waehlngtoi  25»  D.C. 

17*  19  CoBsandlng  General,  Continental  Army  Connand,  Ft. 
Monroe,  Va. 

20  President,  Board  #1,  Headquarters,  Continental  Army 

Command,  Ft.  Sill,  Okla. 

21  President,  Board  #2,  Headquarters,  Continental  Army 

Command,  Ft.  IQiox,  Ky. 

22  President,  Board  #3,  Headquarters,  Continental  Army 

Coornmnd,  Ft.  Banning,  Ga. 

23  President,  Board  #4,  Headquarters,  Continental  Army 

Connand,  Ft.  Bliss,  Tex. 

24  CooBandlng  General,  U.S.  Army  Caribbean,  Ft.  Amador, 

C.'w.  ATTN:  Cnl.  Off. 

25  Cosmnrder-ln-Chlef , European  Coamtand,  APO  126,  c/o  FM, 

Kev  7ork,  N.T. 

26-  27  CooBander-ln-Chlef , Far  East  Conand,  APO  500,  c/o  B4, 
San  Francisco,  Calif.  A!^:  ACofS,  J-3 
26-  29  Comandlng  General,  C7.S.  Army  Europe,  APO  403,  c/o  Bi, 
New  Tork,  N.T.  ATOf:  OFOT  Dlv.,  Combat  Dev.  Br. 

30-  31  Coasandant,  Conand  and  General  Staff  College,  Ft. 
Teavaiivorth,  Kan.  ATOI:  A1XI£(AS) 

32  CoiusridSLnt,  The  Artillery  and  Guided  Missile  School, 

..  Gill,  Oklh. 

33  •’Cs.ratary,  The  Antiaircraft  Artillery  and  Guided 

Missile  School,  Ft.  Bliss,  Texas.  ATTN:  Maj. 

George  L.  Alexander,  Dept,  of  Tactics  and 
Combined  Arms 

34  Coaendlng  General,  Medical  Field  Service  School, 

Brooke  Army  Medical  Center,  Ft.  Sam  Houston,  Tex. 

35  Director,  Special  Weapons  Development  Office, 

Headquarters,  CONAPC,  Ft.  Bliss,  Tex.  ATTN:  Lt. 

Arthur  Jasklerny 

36  CoMandant,  Army  Medical  Service  Graduate  School, 

Walter  Reed  Army  Medical  Center,  Washington  25i  D.C. 

37  Superintendent,  U.S.  Military  Academy,  West  Point,  N.T. 

ATTN:  Prof,  of  Ordnance 

38  ri—Miiiilaiit . Chemical  Corps  School,  Chemical  Corps 

Training  Coaand,  Ft.  McClellan,  Ala. 

39  CooDsndlng  General,  Besesreh  and  Engineering  Coi^nd, 

Army  Chemical  Center,  Md.  ATTN:  Deputy  for  RW  and 
Non-Toxic  Material 

40-  41  Cooendlng  General,  Aberdeen  Proving  Grounds.  Md. 

(inner  envelope)  ATOI:  RD  Control  Officer  (for 
Director,  Ballistics  Research  Laboratory) 

42-  44  Conandlng  General,  The  Engineer  Center,  Ft.  Belvolr, 
Ta.  AtnE:  Aset.  Coemandant,  Engineer  School 

45  CoHMndlng  Officer,  Engineer  Research  end  Development 

Laboratory,  Ft.  Belvolr,  Va.  AnSf:  Chief,  Technical 
Intelligence  Branch 

46  Ccamandlng  Officer,  Plcatlnny  Arsenal,  Dover,  N.J. 

Ansi:  ORDBB-TR  ' 


47  CoB&anding  Officer,  Anny  Medical  Research  Laboratory, 
Ft.  Knox,  Ky. 

46-  49  Commanding  Officer,  Chemical  Corps  Chemical  and  Rudlo- 
logical  Laboratory,  Army  Chemical  Center,  Md.  ATTN: 
Tech.  Library 

50  Cocmaandlng  Officer,  Transportation  RIcD  Station,  Ft. 

Bustle,  Va. 

51  Director,  Technical  Documente  Center,  Evans  Signal 

Laboratory,  Belmar,  N.J. 

52  Director,  Watervaye  Experiment  Station,  PO  Box  631, 

Vicksburg,  Miss.  ATTN:  Library 

53  Director,  Armed  Forces  Institute  of  Pathology,  7th  and 

Independence  Avenue,  S.W.,  Washington  2%  D.C. 

54  Director,  Operations  Research  Office,  Johns  Hopkins 

University,  7100  Connecticut  Ava.,  Chevy  Chase,  Md. 
Washington  15,  D.C. 

55*  56  Conmandlng  General,  Quartermaster  Research  and 
Development  Command,  Quartermaster  Research 
and  Development  Center  Natick,  Mesa.  ATIN: 

CBR  Liaison  Officer 

57-  63  Technical  Information  Service,  Oak  Rtdge,  Tenn. 
(Surplus) 


NAVT  ACTIVITIES 

64-  65  Chief  of  Naval  Operations,  D^,  Washington  25,  D.C 
ATTN:  OP-36 

66  Chief  of  Naval  Cperatloos,  D/ffg  Waablagton  25,  D.C. 

ATTN:  0P-03BG 

67  Director  of  Naval  Intelligence,  D/N,  Washington  25. 

D.C.  AnN:  0P-922V 

66  Chief,  Bureau  of  Medicine  and  Surgery,  D/N,  Washlngtcr 
25,  D.C.  A^N:  Special  Weapons  Defense  Dlv. 

69  Chief,  Bureau  of  Ordnance,  D/^,  Washington  25,  D.C. 

70  Chlsf,  Bureau  of  Ships,  D/V,  Wae.hL-.gtor  25,  D.C.  ATIN; 

Code  346 

71  Chief,  Biu'eau  of  Yards  and  Docks,  D/N,  Wasliltigtcr.  25. 

D.C.  ATTN:  D-440 

72  Chief,  Bureau  of  Supplies  and  Accounts,  D/N,  Washing- 

ton 25,  D.C. 

73-  74  Chief,  Bureau  of  Aerw..-u*lcs,  D/N,  Wsshl'-.gtoii  25,  D.C. 

75  Chief  cf  Naval  Research,  De'^artaant  of  tne  Navy 

Washington  25,  D.C.  ATTN:  Code  «11 

76  Commander-In-Chief,  U.S.  Pacific  Fleet,  Fleet  Post 

Office,  San  Francisco,  Calif. 

77  CG«ander-ln-Chlef,  U.S.  Atlantic  Fleet,  U.S.  Naval 

Base,  Norfolk  11,  Va. 

78-  61  CoBm»ndQnt,  U.S.  M'lrlne  Corps,  Washlngtor.  25.  D.C. 

ATTN:  Code  A03H 

President,  U.S.  Naval  War  College,  Newport,  R... 

33  Superintendent,  U.S.  Naval  Foett7aduate  School, 
Monterey,  Calif. 

84  CoBandlng  Officer,  U.S.  Naval  Schools  Conand,  U.S. 
Naval  Station,  Treasure  Island,  San  Francisco, 

Calif, 

65  CwBBandlng  Officer,  U.S.  Fleet  Training  Center,  Naval 
Base,  Norfolk  11,  Va.  ^ITN:  Special  Weapone  School 
06-  87  Comandlng  Officer,  U.S.  Fleet  Training  Center,  Naval 
Station,  San  Diego  36,  Calif.  ATTN;  (SWP  School) 
Coonandlng  Officer,  Air  Devalopoent  Squadron  5,  VK-^, 
U.S.  Naval  Air  Stotlon,  Moffett  Field,  Ca\lf, 

89  Conmandtng  Officer,  U.S.  Novul  Damage  Control  Training 

Center,  Navol  Base,  Phlladelphlh  12,  Pa.  ATTN:  ABC 
Defenae  Course 

90  Conaanding  Officer,  U.S.  Nnval  Unit,  Chemical  Corps 

School,  Army  Chemical  Training  Center,  Ft.  McClell-.n, 
Ala. 

91  Coimnander,  U.S,  Naval  Ordnance  Laboratory,  Silver 

Spring  19,  Md.  ATTN:  EB 


UNGLASSIHED 


92  Coonandor,  U.S.  Naval  Ordnance  Laboratory,  Silver 

Spring  19,  Md.  ATTN:  EH 

93  Comander,  U.S.  Naval  Ordnance  Laboratory,  Silver 

Spring  19,  Md.  ATIN:  R 

9^  Conander,  U.S.  Naval  Ordimnce  Teat  Station,  Inyokern, 
China  Lake,  Calif. 

9^  Offlcer-ln^Charge,  U.S.  Naval  Civil  Engineering  Rea. 
and  Evaluation  Lab.,  U.S.  Naval  Conatructlon  Bat* 
tallon  Center,  Port  Bueneae,  Calif.  ATTN:  Code  7^3 

96  CoBBandlng  Officer,  U.S.  Naval  Medical  Research  Inst., 

National  Naval  Medical  Center,  Bethesde  Ik,  Md. 

97  Director,  Naval  Air  Experloental  Station,  Air 

Material  Center,  U.S.  Naval  Baee,  nilladelj^la, 

Penn. 

96  Director,  U.S.  Naval  Research  Laboratory,  Waehlngton 
85»  D.C.  ATTN:  Code  2029 

99  CoOBandlng  Officer  and  Director,  U.S.  Navy  Ilectronlce 
Laboratory,  San  Diego  ^2,  Calif.  ATTN:  Code  k223 
100-101  Cocnanding  Officer,  U.S.  Navel  Radiological  Defenee 

Laboratory,  San  Francisco  2k,  Calif.  ATTN:  Technical 
Inforaatlon  Division 

102-103  CoHundlng  Officer  and  Director,  David  V.  Taylor  Model 
Basin,  Waehlngton  7,  D.C.  kTOt:  Library 
10k  CoMender,  U.S.  Naval  Air  Developaent  Center,  Johne- 
vllle.  Pa. 

10^  Director,  Office  of  Naval  Research  Branch  Office,  1000 
Geary  St.,  San  Pranclsco,  Calif. 

106-111  Technical  Information  Service,  Oak  Ridge,  Tenn. 
(Surplua) 


AIR  FORCE  ACTIVITIES 

112  Asst,  for  Atonic  Energy,  Headquarters,  USA/,  Waehlng- 

ton 25,  D.C.  ATTN;  DCS/O 

113  Director  of  Operations,  Headquarters,  USAF,  Waehlngton 

25,  D.C.  A'nN:  Operations  Analysis 
Ilk  Director  of  Plans,  Headquarters,  USAF,  Washington  25, 
D.C.  ATTN;  War  Plana  Dlv. 

115  Director  of  Research  and  Developaent,  Headquarters, 
USAF,  Washington  25,  D.C.  ATTN:  Coabat  Coaponents 
Dlv. 

116-117  Director  of  Intelligence,  Headquarters,  USAF,  Washing- 
ton 25,  D.C.  ATTN;  AF0rN-IB2 

lid  The  Surgeon  G«neralf  Headquarters,  USAF,  Vashiiigton  25/ 
D.C.  ATTN;  Bio.  Def.  Br.,  Pre.  Med.  Dlv. 

119  Deputy  Chief  of  Staff,  Intelligence,  Headquaz^ters,  U.S. 

Air  Forces  Europe,  APO  633,  c/o  HI,  Nev  Tork,  N.7. 
A^N:  Directorate  of  Air  Targets 

120  Comander,  k9Tbh  Reconnaissance  Technical  Squadron 

(Au0Bented),  APO  633,  c/o  m,  Nev  Tork,  N.T. 

121  CoBBiander,  Far  East  Air  Forces,  APO  925,  c/o  FM,  San 

Francisco,  Calif. 

122  Coanender-in-Chlef , Strategic  Air  Coanand,  Offutt  Air 

Force  Baee,  Omaha,  Nebraska.  ATTN:  Special  Weapons 
Branch,  Inspector  Dlv.,  Inspector  General 

123  Cooaandar,  Tactical  Air  Conaand,  Langley  AFB,  Va. 

Aral:  Doctjaents  Security  Branch 
12k  Connander,  Air  Defenee  Comnand,  Ent  AFB,  Colo. 

125-126  CoLSBander,  Wright  Air  D«velof«ent  Center,  Wrlght- 
Pattaraon  AFB,  Dayton,  0.  ATTN;  WCRSN,  Blast 
Effects  Research 

127  Comander,  Air  Training  Ccmaand,  Scott  AFB,  Belleville, 

111.  ATTN:  DCS/0  GTP 

128  Assistant  Chief  of  Staff,  Installations,  Headquarters, 

USAF,  Washington  25,  D.C.  ATTN:  AFCIE-E 

129  Comandar,  Air  Research  and  Developtoent  Cosmand,  PO 

Bo*  1395,  Baltimore,  Md.  ATTN:  RDDN 

130  CoOBander,  Air  Proving  Ground  Coonand,  Eglln  AFB, 

Fla.  ATTN:  AG/TRB 

131-132  Director,  Air  University  Library,  Msxwell  AFB,  Ala. 
133"!*^  Coanander,  Flying  Training  Air  Force,  Waco,  Te*. 

ATTN;  Director  of  Observer  Training 
Ikl  CoimBander,  Crew  Training  Air  Force,  Randolph  Field, 

Te*.  ATTN;  2CTS,  DCS/o 


lk2  Comnder,  Headquarters,  Technical  Training  Air  Force, 
Gulfport,  Mies.  ATTN:  TA&D 

lk3-lkk  Conaandunt,  Air  Force  School  of  Aviation  Medicine, 
Randolph  AFB,  Tex. 

Ik5*150  CoQBaiuier,  Wrlglit  Air  Developosnt  Center,  Wrlgbt- 
Patterson  AFB,  Dayton,  0.  ATTN:  WCOSI 
151-152  Coonander,  Air  Force  Cambridge  Research  Center,  ISj 
Hanacom  Field,  Bedford,  Mast.  Attn:  C1^3T-2 
X53*'155  CooBander,  Air  Force  Special  Weapons  Center,  Klrtland 
AFB,  N.  Me*.  ATW:  Library 

156  Coom'indant,  USAF  Institute  of  Teclmology,  Wrlght- 

Pattersoii  AFB,  Dayton,  0.  ATTN:  Resident  College 

157  CoBBH.ndar,  Lowry  AFB,  Dat.ver,  Colo.  ATTN:  Department 

of  Armament  Tr*ilnli<g 

158  Commander,  100yu<  Special  Weapons  Squadron,  Head- 

quarters, USAF,  W'lshlngton  25,  D.C. 

159-160  INe  RAND  Corporation,  1700  Main  Street,  Santa  Monica, 
Cellf.  ATHI:  Nuclear  kr.er^  Division 

161  Conoander,  Second  Air  Force,  Barksdale  AFB,  Louisiana. 

ATTN:  Operations  Analysis  Office 

162  Connander,  Clgt.th  Air  Force,  Westover  AFB,  Maes  ATTN; 

Operations  Analysis  Office 

163  fOBnander,  Fifteer.th  Air  Force,  March  AFB,  Calif. 

ATTN:  Operations  Analysis  Office 
l6k-IT0  Technical  I{.l'0'matlor.  r>ervlce,  Oak  Ridge,  Tenn. 
(Surplus) 


OTHER  DEPARTKERT  Of  DEFENSE  ACTXVmBS 

171  Aast.  Secretary  of  Defense,  Research  and  Development, 
d/d,  Washington  25,  D.C.  ATTN:  Tech.  Library 

1T2  U.S.  Doc..nBnte  Officer,  Office  of  the  U.S.  National 
Military  Representative,  SHAPE,  APO  55,  Naw  York, 

N.Y. 

173  Director,  Weapons  Systems  Ev&l\iatlon  Group,  OSD,  Rm 
2E1006,  Pentagon,  Washington  25,  D.C. 

17k  Armed  Services  £*plOBlvea  Safety  Board,  D/D,  Building 
T-7,  Gravelly  Point,  Washington  25,  D.C. 

175  CoMandant,  Armed  Forces  Staff  College,  Norfolk  11, 

Ve.  Aral:  Secretary 

176-ldl  rifiniling  Genmral,  Field  CaMaDd,  Armed  Forcee  Spe- 
cial Weapons  Project,  FO  Box  5100,  Albuquerque,  I. 

Me*. 

182-183  CoMuidlng  General,  Field  Ccnaad,  Armed  Forces,  Special 
Weapons  Project,  PO  Bo*  5100,  Albuquerque,  N.  Hex. 
Aral:  Technical  Training  Grotq> 

l6k-192  Chief,  Armed  Forces  Special  Weapons  Project,  Washington 
25,  D.C.  ATTN:  Documents  Library  Branch 

193  Office  of  the  Technical  Director,  Directorate  of  Ef- 
fects Tests,  Field  Ciwand,  AFSWP,  FO  Box  577, 

Menlo  Park,  Calif.  AnS:  Dr.  E.  B.  Doll 

194.200  Technical  Information  Service,  Oak  Ridge,  Tenn. 

(Surplus) 


ATOMIC  BIERCr  CQMtlSSIOR  ACTIVITIES 

201-203  U.S,  Atomic  Energy  Comlaelon,  Claaalfled  Technical 

Library,  I9OI  Constitution  Ave.,  Washington  25,  D.C. 
Aral:  Mrs.  J.  M.  O'Laary  (For  WA) 

20k-205  Los  Alaatos  Scientific  Laboratory,  Report  Library,  PO 
Box  1663,  Los  Alamos,  N.  Mex.  Aral:  Helen  Redman 
206-210  Sandta  Corporation,  Classified  Document  Division, 
Sandla  Base,  Albuquerque,  N.  Mex.  Aral:  Martin 
Lucero 

211-215  University  of  California  Radiation  I.aboratory,  PO  Bo* 
806,  Livermore,  Calif.  ATTN:  Margaret  Edlund 
21k  Weapon  Data  Section,  Technical  Information  Service, 
Oak  Ridge,  Tenn. 

215-274  Technical  Inforaatlon  Service,  Oak  Ridge,  Tenn. 

( Surplus ) 


